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ABSTRACT 
The  extreme  speed  and  altitude  ranges  of  which  modern,  high  perfor- 
mance, supersonic  airplanes  are  capable  have  become  manifested  in  wide 
variations  in  airframe  transfer  function  parameters  as  flight  conditions 
change.   Airframe  control  systems  and  stability  augmentation  devices 
utilized  in  these  aircraft  require  some  form  of  gain  changing  to  pro- 
vide satisfactory  stability  and  control  characteristisc  over  the  complete 
flight  envelope.   In  present  day  operational  aircraft,  system  gains  are 
scheduled  as  functions  of  airspeed  and  altitude,  however  these  mechaniza- 
tions have  seme  shortcomings  which  are  pointed  out  in  this  paper,  and  also, 
the  desirability  of  self-adaptive  aircraft  control  systems  is  discussed. 
A  survey  of  adaptive  techniques  which  have  been  applied  to  the  airplane 
control  problem  is  presented,  and  a  self-adaptive  pitch-stability  aug- 
mentor  which  utilizes  a  digital  computer  in  the  adaptation  loop  is  pro- 
posed.  The  proposed  self-adaptive  system  is  tested  in  conjunction  with 
a  digital  three  degree  of  freedom  flight  simulation  of  the  pitch  dynamics 
of  a  representative  supersonic  fighter-attack  aircraft,  and  is  found  to  be 
capable  of  performing  the  adaptation  process  in  two  cycles  of  airplane 
short-period  oscillation.   Requirements  for  practical  mechanization  in 
a  flight  vehicle  and  recommendations  for  future  research  are  discussed. 
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T  -  INTRODU 

As  the  speed  and  altitude  capabilities  of  modern  high  performance 
supersonic  aircraft  have  increased  in  recent  years,  the  range  of  varia- 
tion of  the  dynamic  characteristics  of  these  aircraft  have  increased  pro- 
portionately.  A  ten  to  one  variation  in  damping  ratio  and  natural  fre- 
quency of  some  of  the  dominant  modes  of  aircraft  motion  is  not  uncommon 
for  modern  airplanes  with  Mach  2  and  70,000  ft.  altitude  capabilities. 
As  an  example,  for  the  airplane  considered  in  this  study,  the  damping 
ratio  of  the  dominant  longitudinal  mode  varies  f rom  £   -   0.445  at  M  =  0.9 
@  10,000  ft.  to£=  0.064  at  M  =  1.2  (?  45,000  ft. 

Various  stability  augmentation  devices  have  been  used  to  decrease  the 
range  of  variation  of  airplane  parameters  and  to  hold  them  as  near  as 
possible  to  optimum  values.   These  devices,  in  their  historical  order  of 
application,  include  fixed  gain  dampers,  pilot-adjusted-gain  dampers, 
speed-and-altitude-programmed-gain  dampers,  and  self-adaptive-gain-adjust- 
ing dampers.   A  number  of  methods  have  been  investigated  for  mechanizing 
self-adaptive  stability  augmentation  systems  for  aircraft,  utilizing  vari- 
ous techniques  which  will  be  discussed  later.  And,  the  purpose  of  this 
study  is  to  investigate  the  use  of  an  airborne  stored  program  digital 
computer  for  generation  of  the  adaptation  commands  in  a  self-adaptive 
pitch  damping  system  to  be  used  in  a  high  performance  fighter- interceptor 
type  aircraft.   The  dynamics  of  the  airplane,  including  control  system 
components,  and  the  adaptation  computer  were  simulated  on  a  Control  Data 
Corporation  1604  general  purpose  digital  computer. 

Before  beginning  the  discussion  of  adaptive  control  systems  as  applied 
to  aircraft,  it  is  felt  that  a  better  understanding  of  the  system  require- 
ments may  be  gained  by  a  short  review  of  the  characteristic  dynamic  modes 
of  a  supersonic  aircraft  and  the  theory  underlying  the  analysis  and  syn- 
thesis of  these  motions. 


F  THE  i 
The  airplane      axis  system  which  is  generally  used  in  airplane 
stability  and  contr.     --lies  is  shown  in  figure  1,  as  three  mutually 
perpendicular  axes  intersecting  at  the  airplane  center  of  gravity.  The 
longitudinal,  or  X,  axis  lies  in  the  airplane's  vertical  plane  of  sym- 
metry and  is  aibitrarily  oriented  by  the  airplane  designer  so  as  to  de- 
fine a  horizontal  plane  of  symmetry  of  the  fuselage.   In  almost  all  air- 
craft the  X  body  axis  does  not  coincide  with  the  longitudinal  principal 
axis  (longitudinal  axis  about  which  the  moment  of  inertia  is  least)  due 
to  the  high  location  of  the  mass  of  the  tail  surfaces  and  due  to  the 
assymetric  top-to-bottom  configuration  of  the  fuselage.   Consequently,  the 
longitudinal  principal  axis  normally  slopes  downward  from  tail  to  nose  and 
coincides  with  the  X  body  axis  only  at  the  center  of  gravity.  The  dynamic 
unbalance  of  mass  about  the  X  axis  thus  causes  the  airplane  to  have  an 

I  „  product  of  inertia  which  will  be  important  when  cross  coupling  of  air- 
xz 

plane  motions  is  discussed. 

The  Z  body  axis,  or  normal  axis,  is  perpendicular  to  the  X  axis  and  to 
the  wing  plane  and  has  positive  direction  downward.   And,  the  Y,  or  lateral 
axis,  is  simply  defined  as  being  normal  to  the  X  and  Z  axes.   Fortunately, 
due  to  the  lateral  symmetry,  an  airplane  has  no  dynamic  unbalance  about  the 
Y  axis  and  thus  only  one  product  of  inertia,  I   ,  exists. 

A  second  set  of  references,  the  earth  reference  system,  must  also  be 
considered  due  to  action  of  components  of  gravity  along  the  three  axes. 
To  determine  the  gravity  components,  the  airplane's  attitudes  in  pitch  and 
roll  must  be  determined  and  the  acceleration  of  gravity  resolved  along  the 
three  axes  by  utilizing  the  Euler  equations.  The  pitch  and  roll  angles  are 
illustrated  in  figure  2,  where  pitch  angle  is  shown  with  wings  level  and 
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Figure  1.   Body  axis  system  and  parameters  used  in  the 
digital  simulation 


toll  Si  In  these  two  conditions  the  earth 

reference  angles  are  equivalent  to  integrations  of  airplane  angular  vel- 
ocities about  their  respective  body  axes.   But,  if  the  craft  is  in  an 
attitude  with  non-zero  values  of  both  pitch  and  roll,  these  angles  must  be 
obtained  from  the  projection  of  the  X  and  Y  axes  onto  two  vertical  planes 
containing  these  axes. 

A  third  earth  reference  angle  which  has  not  been  illustrated  and  which 
plays  no  part  in  the  dynamics  of  motion  is  the  yaw,  or  heading,  angle  which 
is  the  projection  of  the  longitudinal  body  axis  onto  the  horizontal  plane. 

Characteristic  Airplane  Dynamics 
The  six  degrees  of  freedom  of  motion  of  an  aircraft  can  now  be  defined 
as  three  angular  freedoms  in  roll  pitch  and  yaw  about  the  X,Y,  and  Z  axes 
respectively;  and  longitudinal,  lateral,  and  normal  translations  along  the 
X,  Y,  and  Z  axes  respectively.  Thus,  the  dynamics  can  be  represented  in 
terms  of  six  simultaneous  differential  equations  of  motion.   Before  dis- 
cussing the  makeup  of  the  six  equations  of  motion,  it  is  felt  that  a  better 
understanding  of  the  following  material  may  be  obtained  by  a  review  of  the 
various  modes  of  aircraft  motion.  The  motions  can  generally  be  divided  in- 
to two  classes,  longitudinal  modes  and  lateral-directional  modef.   The  first 
class  involves  pitch  rotation  combined  with  translations  along  the  X  and  Z 
axes,  whereas  the  lateral-directional  modes  involve  rollings  yawing,  and 
translation  along  the  Y  axis.  This  division  of  modes  occurs  because  of  the 
low  level  of  aerodynamic  and  inert ial  cross  coupling  between  longitudinal 
and  lateral  motions.   For  small  amplitude  oscillations  (excursions  in  pitch, 
roll,  and  yaw  less  than  15-20  degrees)  the  two  classes  can  be  assumed  dynami- 
cally independent. 
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Figure  2.  Earth  reference  system  used  in  the  digital  simulation 


i  lateral 
directional  modes         ssentative  supersonic  airplane.   The  top  half 
of  the  S  plane  is  shown  for  each  case.   It  is  seen  that  the  longitudinal 
motions  occur  in  two  oscillatory  modes,  one  with  a  relatively  long  peri' 
the  phugoid  mode,  and  one  with  a  much  shorter  period,  the  short-period 
mode.   The  damping  of  the  phugoid  is  usually  low  and  in  some  cases  it  is 
a  neutral  or  even  divergent  oscillation.  However,  due  to  its  long  period, 
the  motion  is  easily  kept  under  pilot  control  with  little  or  no  conscious 
attempt  at  correction,  and  no  artificial  stability  augmentation  is  re- 
quired for  this  mode.   The  phugoid  occurs  at  constant  angle  of  attack  and 
manifests  itself  as  a  continuous  trade-off  between  airspeed  and  alt  tude. 

Figures  4,  and  5.  show  that  the  period  and  damping  ratio  of  the  short 
period  mode  of  a  supersonic  aircraft  can  vary  over  wide  limits,  and  the  in- 
clusion of  the  time  lag  of  the  human  pilot  in  the  control  loop  can  some- 
times cause  the  short-period  poles  to  move  into  the  right  hand  side  of  the 
S  plane.  This  necessitates  the  use  of  pitch  stability  augmentation  not 
only  for  improvement  of  handling  qualities  but  for  reasons  of  safety. 
This  mode  occurs  at  constant  airspeed  and  appears  as  a  pitching  oscilla- 
tion accompanied  by  small  amplitude  translations  up  and  down  the  Z  axis. 
It  is  extremely  important  in  maneuvering  flight  and  is  of  prime  interest  to 
designers  of  airborne  fire  control  and  stability  augmentation  systems. 

The  lateral-directional  pole-zero  plot  for  a  supersonic  aircraft 
normally  exhibits  only  one  pair  of  complex  poles  and  two  real  poles  as 
shown  in  fig.  3.  The  frequency  and  damping  ratio  of  the  lateral-direction- 
al short  period,  or  Dutch  roll,  mode  vary  as  flight  conditions  change  in  a 
manner  similar  to  the  variation  of  longitudinal  short-period  characteristics. 
This  mode  is  a  combined  rolling -yawing  oscillation  with  s      mplitude 
lateral  translations,  ^.nd  Its  name  is  derived  from  the  swaying  motions  of 
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Figure  3   Pole-Zero  Plot  of  the  Characteristic  Modes  of 
Motion  of  a  Supersonic  Aircraft 
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a  Dutch  skater.   The  Dutch  roll  oscillations  of  a  supersonic  aircraft 
usually  require  artificial  clamping  ir/  the  same  manner  as  required  by  the 
longitudinal  short-period  mode.   Therefore,  even  though  this  study  is  con- 
cerned with  an  adaptive  longitudinal  stability  augmentation  system,  the 
results  and  conclusions  are  applicable  to  the  lateral-directional  case 
as  well. 

The  roll  damping,  or  roll  subsidence,  mode  is  shown  as  a  negative 
real  pole  and  it  represents  the  viscous  roll  damping  of  the  aircraft. 
And,  the  spiral  mode  is  usually  a  steady  divergence  which  manifests  itself 
as  a  gradually  steepening  turn  with  accompanying  sideslip.   The  divergence 
is  fortunately  slew  and,  as  in  the  case  of  the  longitudinal  phugoid  mode, 
is  easily  controlled  by  the  pilot.   The  roll  subsidence  and  spiral  mode 
poles  can  occasionally  form  a  complex  pole  pair  in  the  case  of  some  sub- 
sonic airplanes  with  large  dihedral  effect.   This  lateral-directional 
"phugoid"  mode  will  exhibit  small  amplitude,  long  period  rolling  oscillations 
with  considerable  lateral  translation  and  very  little  yawing  oscillation. 

This  concludes  the  survey  of  normal  modes  of  airplane  motion.   How- 
ever, there  are  other  "patha logical"  modes  which  can  be  experienced  when 
the  rates  and  amplitudes  of  motion  become  large.   Some  of  these  are  the 
phenomena  of  inertial  cross  coupling  in  roll,  autorotation  (Ref.  125.), 
and  longitudinal  pitch-up  (Ref.  86);  all  of  which  can  be  eliminated  or  al- 
leviated by  the  use  of  well  designed  airplane  stability  augmentation  systems. 
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Aerodynamic  Stability  Derivatives 
One  final  topic,  the  concept  of  aerodynamic  stability  derivatives, 
should  be  mentioned  before  proceeding  to  the  equations  of  motion.   The 
forces  and  moments  which  act  on  an  aircraft  are  generally  expressed  in 
terms  of  non-dimensional  coefficients  multiplied  by  appropriate  aircraft 
and  flight  condition  parameters.   As  an  example,  the  equation  for  airplane 
lift  is  shown  below. 

Lift  =  CL  §   V  S 
Where : 

Lift  is  in  pounds 

C   is  the  airplane  lift  coefficient 

p  s  atomspheric  density  in  slugs/ft" 

V  =  airplane  forward  velocity  in  ft/sec 

2 
S  =  airplane  wing  area  in  ft  . 

9      2 
The  expression  -r  V     is  a  measure  of  the  kinetic  energy  of  the  air- 
stream  and  is  termed  the  dynamic  pressure  because  it  also  represents  the 
force  which  would  theoretically  be  exerted  against  a  unit  square  surface 
held  normal  to  the  stream.   Thus,  the  product  of  dynamic  pressure  and  wing 
area  represents  a  hypothetical  aerodynamic  force  from  which  the  actual  air- 
plane lift  can  be  obtained  through  multiplication  by  the  lift  coefficient. 
Moments  acting  about  the  aircraft  can  be  obtained  in  the  same  manner  with 
the  inclusion  of  a  moment  arm,  for  dimension  considerations,  as  is  done  in 

the  following  pitching  moment  equation. 

P  2  — 
Pitching  Moment  =>  C  »  VSc 

°  m  2 

where  the  moment  arm  c  is  the  wing  mean  aerodynamic  chord.   Corres- 
ponding expressions  for  drag  force,  side  force,  rolling  moment,  and  yaw- 
ing moment  coefficients  are  included  in  the  table  of  symbols.   It  is  seen 
that  there  exists  an  aerodynamic  coefficient  for  each  of  the  six  degrees 
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of  freedom  of  motion.   These  coefficients  and  their  derivatives  are 
usually  determined  empirically  in  wind  tunnel  and  flight  studies  be- 
cause of  the  inaccuracy  of  analytical  determinations. 

Each  coefficient  can  be  expressed  in  a  Taylor's  series  such  as: 

The  effects  of  the  second  order  and  higher  terms  are  insignificant  in  aero- 
dynamics and  are  dropped.   So,  for  practical  purposes  the  coefficient  can 
be  expressed:     r       (^£l\  _l_  (dQk\  ±-  dCu\         , 

The  partial  derivative  terms  such  as  l-^i}  are  called  aerodynamic  stability 
derivatives  and  are  usually  written  for  brevity  in  the  form  C^cc 

.*.  CL=  CLa.CC   +■  Cto;  <X  +  C/.?  3 
Figure  6,  a  curve  of  lift  coefficient  vs.  angle  of  attack,  can  be  used  to 

illustrate  the  physical  meaning  of  stability  derivatives.   Here, I ^  JZ  J 

is  the  slope,  or  derivative,  of  the  airplane  lift  curve.   It  is  seen  that 

the  curve  is  linear  up  to  <X  &  \2° ,      This  indicates  that  AC  -  C      CC     is 

L    '-cc 

physically  valid  for  small  variations  in  CC  .   The  fact  that  almost  all 
aerodynamic  stability  derivatives  remain  nearly  linear  over  a  similar  range 
of  their  corresponding  variables  enables  the  small  amplitude  motions  of  the 
aircraft  to  be  expressed  as  a  set  of  simultaneous  linear  differential  equa- 
tions. 

How  many  stability  derivatives  must  be  considered  in  the  equations  of 
motion9  Given  the  six  force  and  moment  coefficients  C   ,  C  ,  C   ,  C-.  ,  C  , 


*The  lift  coefficient  C  and  the  drag  coefficient  C  act  perpendi- 
cular and  parallel,  respectively,  to  the  relative  wind,  but  it  can  be 
assumed  with  little  error  that  they  are  also  normal  and  parallel  to  the 

X  body  axis.   Then  C_  =  CT  and  C, ,  =  CL. 

L  D 
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10 
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Angle  of  Attack,  deg. 
Figure  6.   Plot  of  Airplane  Lift  Coefficient  vs.  Angle  of  Attack. 

and  C  and  the  fact  that  each  coefficient  can  be "influenced  by: 

1.  angular  displacement  about  each  axis, 

2.  angular  rate  about  each  axis, 

3.  velocity  along  each  axis,  and 

4.  acceleration  along  each  axis 

there  are  12  possible  derivatives  for  each  coefficient  or  72  in  all.  Also, 
control  surface  deflections  can  add  36  more  terms.  ^Fortunately,  airplane 
symmetry,  which  causes  a  lack  of  cross  coupling  between  longitudinal  and 
lateral  -  directional  modes,  greatly  reduces  the  number.  The  magnitudes  of 
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some  of  the  remaining  derivatives  are  so  small  as  to  make  them  insigni- 
ficant.  As  a  result  there  are  approximately  25  derivatives  considered 
important  enough  to  be  included  in  a  simulation  program.  These  are  sum- 
marized in  Appendix  III.   A  list  of  stability  derivatives  for  a  representa- 
tive supersonic  aircraft  is  tabulated  in  Appendix  I.   The  aircraft  is  the 
same  one  for  which  the  longitudinal  short  period  mode  characteristics  are 
illustrated  in  figures  4,  and  5. 

The  Airplane  equations  of  Motion  and  the  Longitudinal  Short  Period  Transfer 
Function. 

The  six  equations  of  motion  corresponding  to  the  six  degrees  of  free- 
dom of  an  aircraft  are  derived  in  appendix  IV  and  the  final  results  are 
presented  in  table  I.   It  should  be  noted  that  the  equations  represent  the 
balance  of  inertial  and  aerodynamic  forces  and  moments  along  and  about  the 
airplane  body  axes,  and  that  the  terms  are  composed  of  airplane  physical 
constants,  aerodynamic  stability  derivatives,  and  aerodynamic  variables. 
Assumptions  used  in  the  derivation  of  the  equations  are: 

1.  The  change  in  aerodynamic  forces  and  moments  due  to  deviations 
from  an  equilibrium  condition  are  assumed  directly  proportion- 
al to  the  deviation. 

2.  The  deviations  of  the  airplane  from  equilibrium  are  assumed 
small, 

3.  The  airplane  is  a  rigid  body. 

4.  Constant  thrust  is  assumed  during  an  airplane  distrubance, 
and  secondary  power  effects  are  neglected. 

The  equations  are  used  directly  in  the  digital  computer  simulation  of  the 
airplane  dynamics,  which  is  outlined  in  appendix  VI,  and  the  airplane 
longitudinal  short-period  transfer  function  is  derived  from  the  pitching 
and  normal  translation  equations  of  motion  in  appendix  V.   This  transfer 

14 


function  is  of  prime  importance  in  the  study  of  the  airplane  pitch  damper 
because  the  short-period  mode  is  the  dominant  longitudinal  mode.   The 
short-period  transfer  function  is  written  below  both  in  its  general  form  and 
also  in  terms  of  the  airplane  physical  characteristics,  airplane  stability 
derivatives,  and  aerodynamic  parameters.   The  second  form  of  the  transfer 
function  is  shown  to  illustrate  the  dependence  of  the  dynamics  of  the  short 
period  mode  on  changes  in  stability  derivatives  and  aerodynamic  parameters, 
and  it  is  possible  from  inspection  of  the  function  to  determine  which  de- 
rivatives and  parameters  exert  influences  on  a  given  characteristic,  such 
as  damping  ratio,  of   the  short -period  mode. 


Longitudinal  Short  Period  Transfer 
Function 


A.  General  Form 

B.  Detailed  Form 
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Ill  -  T       )F  A  PITCH  DA3        LOMGIT  NATION 

The  purpose  o£  airplane  stability  augmentation  is  to  reduce  the 
range  of  variation  of  a  dynamic  characteristic,  or  characteristics ,  and 
to  hold  the  variables  of  interest  as  close  as  possible  to  predetermined 
optimum  values.   In  the  case  of  an  aircraft  these  optimum  values  are  de- 
fined in  terms  of  "airplane  handling  qualities"  which  describe  an  air- 
plane's response  to  manual  control  by  a  pilot.   Even  though  the  determina- 
tion of  desirable  handling  qualities  is  subjective,  and  is  dependent  on 
pilot  opinion,  a  number  of  studies  have  been  conducted  to  specify  these 
qualities  in  quantitative  engineering  terras.  One  of  the  better  known 
studies,  conducted  at  the  Cornell  Aeronautical  Laboratory,  (Ref,  34),  con- 
cludes that  the  desirable  natural  frequency  and  damping  ratio  of  the  short- 
period  mode  of  a  fighter  airplane  lie  in  the  range  C0n   s  .4  to  .5  cycles  per 
second  and  damping  ratio  =  .65  to  .70.  The  same  study  also  concludes  that 
even  though  the  longitudinal  natural  frequency  is  of  importance  in  maneuver- 
ing flight,  it  can  vary  over  wide  limits  and  the  handling  qualities  will 
still  be  considered  tolerable,  as  long  as  the  short  period  damping  ratio  is 
held  in  the  .65  to  .70  range.   Thus  it  is  considered  more  desirable  to  equip 
a  high  performance  airplane  with  artificial  damping  than  to  provide  static 
stability  augmentation.   Tt  should  be  noted  that  for  the  craft  considered 
in  this  study,  as  well  as  for  almost  all  other  airplanes  of  similar  type 
and  mission,  the  longitudinal  damping  ratio  of  the  un-augmented  airplane 
is  always  below  the  desired  f&   .7,  as  is  shown  in  figure  4.  This  means 
that  the  pitch  damping  system  will  work  to  supplement  the  basic  aerodynamic 
clamping  in  all  cases. 
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Mechanization  of  the  Pitch  Damper 
In  the  discussion  of  the  airplane  longitudinal  short-period  transfer 
function  it  was  shown  that  the  short -period  motion  can  be  represented  as  a 
second  order  system  with  characteristic  damping  ratio  and  natural  frequency. 
It  was  also  shown  that  the  natural  frequency  and  damping  ratio  can  be  ex- 
pressed as  functions  of  the  aerodynamic  stability  derivatives.   The  ex- 
pression for  short  period  damping  ratio  being: 
A     r       _  AC* r 


vfMf  r  A    r         Ac   r„. 

From  the  above  relationship  it  can  be  seen  that  the  most  effective  means 
of  varying  the  damping  ratio  is  to  artificially  vary  the  pitch  damping  de- 
rivative, ^-/mq+oL        *    wn*ch  represents  the  airplane's  aerodynamic  viscous 
resistance  to  pitching  velocity. 

The  pitch  damper  must  then  provide  an  aerodynamic  pitching  moment  which 
is  proportional  to  airplane  pitch  rate,  and  which  has  a  direction  in  opposi- 
tion to  the  angular  rate.  Thus,  the  moment  produced  is  in  phase  with  the 
viscous  aerodynamic  damping  in  pitch  and  the  two  moments  combine  to  form 
the  total  pitch  damping  moment.   The  artificial  damping  moment  is  usually 
generated  by  deflecting  the  airplane's  horizontal  stabilizer  in  direct 
proportion  to  measured  pitch  rate. 

Figure  7.  is  a  block  diagram  of  the  pitch  dampen  ussd  in  this  study. 
The  airplane  transfer  function  shown  in  the  diagram  represents  the  lcn£ 
tudinal  short -period  mode  only  because,  as  will  be  demonstrated  later,  the 
phugoid  mode  and  the  lateral -directional  dynamics  are  of  negligible  impor- 
tance in  longitudinal  stability  augmentation  considerations. 
The  stabilizer  actuator  in  the  diagram  is  an  irreversible  hydraulic  position 
servo,  and  the  parameters  in  the  transfer  function  represent  the  state~o£- 
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Figure  7.   Block  Diagram  of  the  Pitch  Damper 
the-art  as  used  in  modern  aircraft.   It  is  generally  assumed  that  the 
dynamics  of  the  rate  gyro  have  negligible  effect  on  system  operation  be- 
cause the  ratio  of  the  characteristic  frequency  of  state-of-the-art  rate 
gyros  to  the  short-period  mode  frequency  is  of  the  order  of 
(Ref.  43),  and  thus  the  gyro  transfer  function  is  assumed  to  be  unity. 
Table  II  lists  the  11  airplane  transfer  functions  used  in  the  study  along 
with  corresponding  flight  conditions.   And,  figures  8  thru  18  are  roof 
locus  plots  showing  the  behavior  of  the  system  with  varying  damper  gain  at 
each  condition.   The  square  symbols  on  the  plots  indicate  the  gain  values 
at  which  the  airplane  short-period  damping  ratio  has  been  increased  to  0,7. 
It  is  seen  that  the  use  of  the  damper  has  vary  little  influence  on  the 
short  period  frequency  when  operated  over  the  range  of  gains  required  to 
bring  £  up  to  0.7  at  each  flight  condition,   And,  it  is  also  seen  that 
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the  actuator  pole  moves  in   such  a  direction,  with  increasing  damper  gain, 
that  the  actuator  damping  ratio  regains  close  to  0.7  over  the  damper  opera- 
ting range.  The  actuator  mode  will  become  unstable  at  extremely  high  gains, 
but  operation  at  such  gains*  is  of  no  interest  in  this  study. 

It  is  known  that  the  closing  of  a  feedback  loop  around  a  plant  having 
varying  characteristics  will  cause  the  plant  output  to  vary  over  a  smaller 
range,  as  plant  characteristics  vary,  than  the  range  of  variation  of  out- 
put with  no  feedback.   From  this  it  would  be  expected  that  the  pitch  damper 
should  not  only  increase  the  average  value  of  damping  ratio  over  all  flight 
conditions,  but  also  that  the  range  of  variation  of  damping  ratio,  as  condi- 
tions vary,  should  decrease.   This  is  borne  out  in  Figure  19.   which  is  a 
plot  of  shon-period  daraping  ratio  vs.  speed  and  altitude  with  constant 
damper  gain  =  0.3.   The  chosen  value  of  damper  gain  is  approximately  the 
median  of  the  gains  required  for  0,7  damping  ratio  over  all  11  flight 
conditions.   In  the  case  of  the  basic  airplane,  6,       varies  from  0.445  at  M 
=  0.9  @  10,000  ft.  to  £  a  0.064  at  M  =  1.2@  45,000  ft.  whereas  the  varia- 
tion of  /  with  damper  gain  =  0.3  is  from  1.1  to  0.436  for  the  same  two 
flight  conditions.  Thus  the  range  of  variation  about  the  mean  has  been 
reduced  from  approximately  190%  to  approximately  95%. 

Even  though  the  use  of  the  constant  gain  damper  has  caused  a  consider- 
able improvement  in  the  dynamics  of  the  airplane,  a  further  restriction  in 
the  range  of  variation  of  damping  ratio  is  usually  desired  for  fighter-inter- 
ceptor type  aircraft.   One  method  of  solving  this  problem^  which  is  in 
general  use  today,  is  the  programming  of  stability  augmentation  system 
gains,  as  functions  of  measured  flight  parameters,  in  such  a  manner  that 
the  programmed  gains  will  cause  near  optimum  system  response  over  all 
flight  conditions.   The  techniques  of  gain  scheduling  methods  range  from 
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variation  of  damper  gain  as  a  function  of  measured  airspeed  to  systems 
which  program  as  functions  of  airspeed,  altitude,  fuel  load,  center  of 
gravity  position,  wing  flap  position,  landing  gear  posit ion,  and  vari- 
ous combinations  thereof.   One  fact  becomes  obvious  from  the  preceeding 
discussion;  the  dynamic  characteristics  of  the  aircraft  must  be  accur- 
ately known  in  advance  at  all  flight  conditions  and  for  all  airplane 
configurations  if  the  programming  system  is  to  operate  successfully. 
Such  a  detailed  determination  of  airplane  characteristics  requires  a 
large  amount  of  wind  tunnel  and  flight  testing  which  is  both  time  consuming 
and  costly.   And,  even  after  the  desired  information  has  been  determined 
it  is  usually  necessary  to  conduct  lengthy  flight  tests  with  the  program- 
ming and  stability  augmentation  systems  installed  in  the  airplane  for  final 
"tailoring"  of  the  aircraft-system  combination.   In  addition,  accurate  in- 
flight air  data  and  airplane  configuration  information  must  be  continusouly 
fed  into  the  system  during  flight  to  provide  the  basis  for  gain  scheduling. 
The  instrumentation  and  transducers  required  are  generally  costly,  and  are 
usually  complex  enough  that  they  create  reliability  and  fail-safety  problems, 

A  possible  solution  to  these  control  problems  is  to  design  a  system 
which  is  intelligent  enough  to  determine  when  the  airplane's  dynamic  per- 
formance has  changed,  and  in  what  manner  it  has  changed,  and  to  use  this 
information  to  adjust  the  stability  augmentation  system  gains  to  obtain 
optimum  performance  over  all  flight  conditions.  This  is  the  concept  of 
the  self-adaptive  control  system. 
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IV  -  SELF-ADAPTIVE  CONTROL  SYSTEMS 
The  concept  of  adaptive  control  has  been  applied  to  a  wide  range  of 
control  system  types  with  the  result  that  a  number  of  definitions  of  adap- 
tive control  exist.   However,  the  most  often  used  definitions,  and  also  the 
one  which  appears  to  best  fit  the  aircraft  control  problem  is  as  follows: 
An  adaptive  system  is  a  control  system  configuration  in  which 
the  measurement  of  process  dynamics  or  signal  characteristics 
is  utilized  to  automatically  adjust  the  controller  in  an  attempt 
to  achieve  optimum  operation  at  all  times. 
From  this  definition  it  can  be  seen  that  the  adaptive  process  must  per- 
form two  functions,  the  first  of  which  is  the  "identification"  function, 
or  the  measurement  of  process  dynamics  or  signal  characteristics  to  deter- 
mine present  system  state  and  compare  it  with  the  desired  state.  And,  the 
second  is  the  "command"  function  which  is  automatic  changing  of  the  system 
so  that  it  will  approach  a  desired  state.   In  most  adaptive  control  studies 
it  is  assumed  that  the  variation  of  system  characteristics  takes  place  slow- 
ly with  respect  to  the  frequency  of  the  system  mode  of  interest,  and  thus 
the  identification  process  can  be  allowed  a  relatively  long  measurement  in- 
terval during  which  to  determine  system  state.   However,  in  the  ail  craft 
problem,  it  is  desirable  to  perform  the  adaptation  process  relatively  rapid- 
ly in  comparison  to  the  quantity  being  controlled,  in  this  case  the  short 
period  mode  dynamics.   This  is  because  virtually  instantaneous  changes  can 
take  place  in  the  airframe  transfer  function  due  to  release  of  payload  or 
certain  Mach  number  effects  on  stability  derivatives.   For  this  reason, 
one  of  the  objectives  of  this  study  is  the  investigation  of  an  adaptive 
technique  capable  of  rapid  identification  and  command.   Another  considera- 
tion in  mechanization  of  the  adaptive  process  is  the  generally  accepted 
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requirement  ?  :onsidered  truly  adap- 

tive, the  identification  process  can  net  involve  measurement  of  any 

table  which  is  external  to  the  system.   To  illustrate  the  concepts 
and  requirements  discussed  above,  a  block  diagram  of  Gibson's  ideal 
generalized  adaptive  system  (Ref.  76)  is  shown  in  Figure  2G. 
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Figure  20.   Gibson's  Tdeal  Generalized  Adaptive 

System 

Classification  of  Adaptive  Systems 

Tn  order  to  compare  the  approaches  which  have  been  taken  to  the 

design  of  adaptive  control  systems,  Aseltine,  in  reference  4.  has  seperat- 

ed  them  into  5  classes  which,  in  some  cases,  have   a  certain  amount  of 

overlapping,  but  even  so,  they  form  a  convenient  framework  for  comparison, 

Aseltine's  classes  are: 
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1)  Passive  Adaptation:   Systei     ch  adapt  without  automatic 
changes  in  s;        racteristics,  hut  rather  through  de- 

;ti  for  uniform  operation  over  wine  variations  in  environ- 
ment. 

2)  Input  Signal  Adaptation:   Systems  which  automatically  ad- 
just their  parameters  in  accordance  with  input  signal 
characteristics «, 

3)  Extremum  Adaptation:   Systems  which  adjust  to  obtain  a  maxi- 
mum or  minimum  of  some  system  variable. 

4)  System  Variable  Adaptation:   Systems  in  which  the  command 
process  is  actuated  directly  by  measurement  of  a  system 
variable  or  variables. 

5)  System  Characteristic  Adaptation:   Systems  in  which  the 
identification  process  measures  system  variables  in  order 
to  determine  system  characteristics.   The  system  character- 
istics thus  determined  are  compared  \-/ith  desired  character- 
istics and  the  differences  used  to  generate  command  signals. 

It  has  been  noted  from  a  survey  of  adaptive  systems  which  have  been  investi- 
gated for  possible  use  in  aircraft  control  that,  to  date,  all  such  systems 
fit  into  classes  1,  4,  or  5.   And,  it  has  also  been  noted  that  an  alter- 
nate classification  system  might  be  applied  to  aircraft  adaptive  control 
systems.   The  suggested  classes  are: 

1)   High  Gain  Adaptation:   Systems  in  which  the  characteristics 
of  the  aircraft  are  submerged  by  the  use  of  high  feedback 
and/or  forward  loop  gains.   Some  form  of  idealized  model  is 
normally  used  to  provide  either  system  input  signals  or  out- 
put comparison  signals,  and  the  techniques  usually  involve 
either  nonlinear  processes,  with  attendant  limit  cycles  or 
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operation  of  a  minor  system  mode  near  instability,, 

2)  Impulse  Response  Adaptation:   Systems  in  which  the  identi- 
fication process  computes  a  time  history  of  the  plant  im- 
pulse response.  The  impulse  response  is  then  utilized  by 
the  command  process  to  generate  plant  alteration  signals. 

3)  Parameter  Adaptation:   Systems  in  which  plant  parameters 
are  determined  by  direct  computation,,  The  error  between 
the  desired  value  of  the  parameter  and  its  actual  value  is 
used  as  the  command  signal  for  plant  alteration. 

Survey  of  Adaptive  Techniques  which  have  Been  Applied  to  the  Airplane 

Stability  Augmentation  problem 

1)  High  Gain  Adaptation  Systems 

The  earliest  studies  of  aircraft  adaptive  control  involved  systems 
which  operated  at  high  gains  in  order  to  force  the  variation  of  airframe 
parameters  into  a  small  range.   One  of  the  first  of  these  studies  was 
conducted  by  Rath  of  WA.DC  (Ref.  134)  on  a  method  vrhich  was  not  truly  ad- 
aptive, in  that  parameters  were  not  automatically  variedj,  however,  it 
may  be  classed  as  passive  adaptation,   Rath8s,  system  used  a  parallel  com- 
plementary optimum  response  model  as  shown  in  Figure  21.   The  goal  of 
this  study  was  to  design  a  system  which  would  produce  optimum  airframe 
dynamic  characteristics,  as  airframe  parameters  varied,  with  constant 
system  gain.   A  fourth  order  transfer  function  was  selected  as  the  opti- 
mum model  for  the  pitch  dynamics  and  the  same  input  signal  was  applied  to 
the  airframe  and  to  the  model.  The  difference  between  model  and  plant 
output  was  obtained  and  various  functions  of  the  difference  signal  were 
utilized  to  improve  the  airframe  response.   Both  the  forward  loop  gain 
and  the  shaping  loop  gain  were  made  high  in  order  to  "force"  the  system 
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Figure  21.   Rafch's  System  for  Passive  Adaptation  Using  &   Complementary 

Optimum  Response  Model 

overall  characteristics  to  resemble  the  model,  Feedback  functions  studied 
included  pitch  angle,  or  position,  error  feedback,  integrated  position 
error  feedback,  rate  error  feedback,  and  rate-error-plus=>pos  it  ion-error 
feedback.   It  should  be  noted  that  error,  as  referred  to  here,  means  the 
difference  between  plant  output  and  model  output  rather  than  plant  input - 
output  error.   The  first  two  feedback  functions  produced  overall  instabil- 
ity at  high  gains,  whereas  rate  error  feedback  resulted  in  a   stable  system 
with  infinite  steady  state  position  error.   Rate-error-plus-position-error 
feedback  allowed  very  high  gains  to  be  used  before  the  system  became  un- 
stable, and  while  the  system  was  stable,  the  airframe  output  matched  the 
model  output  acceptably  at  most  flight  conditions  while  using  a  constant 
value  of  feedback  gain.   It  was  also  found  that,  whatever  feedback  shaping 
was  used,  the  basic  airframe  and  hydraulic  actuator  must  form  a  stable 
system  with  the  shaped  feedback  loop  alone.   That  is,  the  fe      ?ard  loop 
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ecu VI  not  correct    tabilities  of  '   shaped  feedback 

combination.   The  optimum  transfer  function  feedforward  in  combination 
with  the  shaped  feedl        he  effect  of  changing  the  input  to  the 
airframe  to  a  value  which  caused  the  airframe  to  have  the  desired  out- 
put, assuming  the  airframe  had  the  desired  control  rates  and  high  gain 
control  power. 

An  approach  similar  to  Rath's  passive  adaptation  method  was  taken  by 
Dandois  of  Convair  who  utilized  the  inverse  of  the  optimum  response  model 
in  the  feedback  loop  in  combination  with  a  high  forward  loop  gain  system 
as  shown  in  Figure  22.   (Ref.  43). 
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Figure  22.   Dandois'  System  Using  the  Inverse  of  the  Optimum  Response  Model 

in  the  Feedback  Loop 

He  demon st rated  that,  with  such  a  feedback  arrangement  j,  the  overall  system 

characteristics  approached  the  characteristics  of  the  optimum  response 

model  as  forward  loop  gain  was  increased.   In  actual  operation,  the 

system  held  airframe  parameters  within  acceptable  limits         Lde 


range  o£  flight  conditions,  and  the  high     2  gain  the  better  the 
performance.   Dowevar,  gain  limits,  dependent  upon  flight  conditions, 
were  reached  at  which  the  system  became  unstable.   The  instability  was 
shown  to  be  due  to  sensor  non linearities  rather  than  a  basic,  theoretical, 
system  instability. 

One  of  the  first  "high  gain"  adaptive  systems  to  actually  use  an 
identification  and  command  process  was  the  system-variable-adaptation 
system  of  Flugge-Lotz  and  Taylor  (Ref.  64),  which  is  Illustrated  in  Fig- 
ure 23.   As  is  shown  in  the  block  diagram,  one  of  a  number  of  alternate 
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Figure  23,   Nonlinear  Adaptive  System  of  Flugge-Lot2  and  Taylor 
feedback  paths  was  used  according  to  a  switching  logic  process  based  on 
simultaneous  measurements  of  system  output,  error  signal,  and  derivatives 
of  these.  The  identification  process  thus  consisted  of  the  measurements 
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[back  patbs  and 
the  .".is ted  of  switching  to  the   chosen  path.   Forward  gain  was 

kept  high  in  order  to  "force"  a  small  instantaneous  error  between  input 
and  output,  and  the  system  was  nonlinear  during  any  length  of  time  com- 
parable to  the  natural  period  of  the  plant  because  of  the  switching  be- 
tween feedback  loops  during  one  cycle  of  operation.  The  arrangement  was 
shown  to  be  insensitive  to  variations  in  plant  parameters  over  wide  limits, 
however  a  limit  cycle  persisted  during  operation  due  to  the  feedback  non- 
linearity. 

The  Minneapolis-Honeywell  relay  actuated  adaptive  system  (Ref.  165), 
which  is  presently  undergoing  flight  tests,  is  similar  to  the  Flugge-Lotz 
mechanization  but  uses  only  one  feedback  path  with  switching  logic  in  the 
forward  loop.  The  switching  logic  derives  its  input  from  the  error  signal 
between  optimum  response  model  output  and  the  process  feedback  signal  as 
shown  in  Figure  24.  The  switching  logic  actuates  a  relay,  or  on-off  con- 
troller, to  operate  the  process  at  a  high  gain  with  the  objective  of 
maintaining  a  small  error  signal  between  model  output  and  plaint  output. 
As  in  the  Flugge-Lotz  system,  the  airframe  output  approximately  dupli- 
cates the  model  output  when  the  forward  gain  is  high,  and  the  character- 
istic limit  cycle  is  also  present.  A  second  adaptive  feature  of  this 
system,  aside  from  the  switching  logic,  is  the  use  of  limit  cycle  ampli- 
tude measurement  to  vary  forward  loop  gain.  That  is,  forward  gain  is  ad- 
justed in  such  a  manner  that  the  limit  cycle  amplitude  is  held  constant 
over  all  flight  conditions.  The  prototype  of  this  system  was  successfully 
tested  in  en   F-94C  fighter  aircraft  and  an  improved  version  is  presently 
being  flight  tested  in  the  X-15  high  altitude  research  airplane. 

Another  application  of  the  utilization  of  the  inverse  optimum  re- 
sponse transfer  function  in  the  feedback  loop,  in  conjunction  with 
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Figure  2'+.   Minneapolis -Honeywell  Relay  Actuated  Adaptive  System 
high  forward  loop  gain,  is  the  General  Electric  system  (Ref.  32),  which 
was  recently  tested  in  an  F- 102  interceptor  aircraft.   The  mechanization 
is  similar  to  that  of  L'andois  shown  in  Figure  22,  except  that  an  active 
system-variable-adaptafcion  feature  has  been  added.  This  feature  involves 
variation  of  the  forward  loop  gain,  as  flight  conditions  change,  to 
maintain  the  system  just  below  the  boundai*y  of  instability.   As  was  p 
viously  mentioned,  when  considering  systems  of  this  type,  the  higher  the 
gain,  the  closer  tvill  be  the  airframe  characteristics  to  the  optimum  re- 
sponse characteristics  and  thus  if  the  system  is  operated  near  the  in» 
stability  boundary, the  maximum  performance  will  be  extracted.   In  opera- 
tion, the  adaptation  c,  stem  monitors  the  frequency  of  the  control  surface 
actuator  because  the  actuator  mode  is  the  one  which  becomes  unstable  as 
forward  gain  is  increased.  The  gain  is  then  adjusted  to  r        :  he 
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frequency  of  this  mode  constant  over  all  flight  conditions.   It  should 
be  noted  that  the  actuator  mode  is  a  lightly  damped  oscillation  and 
that  normal  control  inputs  are  utilized  to  provide  the  required  pertur- 
bations. The  frequency  measurement  is  accomplished  with  a  narrow  band 
pass  filter  and  zero  crossing  counter  arrangement. 

2}   Impulse  Response  Adaptation  Systems 

The  function  of  the  identification  process  in  the  impulse  response 
adaptation  class  of  systems  is  first,  the  determination  of  a  time  history 
of  the  plant  impulse  response  through  measurement  of  system  input  and  out- 
put signals  over  a  period  of  time,  and  second,  the  determination  of  plant 
characteristics  through  examination  of  the  impulse  response.   An  early 
method  of  impulse  response  self  optimization  which  was  investigated  by 
Aseltine  (Ref.  4}  serves  to  illustrate  the  process,  even  though  this 
particular  mechanization  would  not  be  suitable  for  aeronautical  use. 
The  impulse  response  was  obtained  by  periodic  input  pulses  which  were 
widely  spaced  with  respect  to  the  period  of  the  plant  mode  of  interest. 
After  each  perturbation,  the  system  damping  ratio  was  computed  by  separ- 
ate integration  of  the  positive  and  negative  portions  of  the  impulse  re- 
sponse and  obtaining  the  ratio  of  positive  to  negative  area.   Damping 
ratio  was  computed  from  this  ratio  and  the  system  damping  coefficient  was 
varied  in  a  closed  loop  manner  to  cause  actual  damping  ratio  to  approach 
a  desired  value.   This  method  was  considered  undesirable  for  aircraft  use 
primarily  because  normal  operating  signals  must  be  excluded  from  the  in- 
put during  the  comparatively  lengthy  identificatior  period. 

A  method  of  impulse  response  determination  in  which  normal  input 
signals  could  be  included  in  the  identification  process  was  investigated 


by  Anderson  .cs  te  t  ion  made      F   the 

convolution  theorem  whereby  che  STrStero  input  convolved  in  time  with  the 
impulse  tespoose  re  time  hist       the  system  output.   In  opera- 

tion, random  noise  w      tocorrelation  functional"; (T)   was  superimposed  on 
the  normal   iput  signal  as  shown  in  Figure  25.  The  plant  output  was  filt- 
ered to  r.-:;  '  ■  |  '^enfc  from  the  total  i     sea  and  the  ii 
put  and  output  noise  components  ^ere  cross-correlated  to        ty)o(T)- 
JjO()(p;((T-X)cj/(      .  The  noise  input  had  a  bandwidth  : 

-CD 

bandwidth  so  that  <pii(T)   was  effectively  an  impulse     (pio(T) 
proximately  the  pi       ilse  response.  Twelve  channels  of  corss  correla- 

n  were  used  and  each  channel  provided  one  point  on  a  time  history  of 
the  impulse  response.  Damping  ratio  was  determined  and  utilized  as  dis- 
cussed in  the  proceeding  paragraph. 
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other  method  oc      se  response  determination  utilizing  the  convolu- 
tion theorem  in  conjunction  with  simultaneous  input  and  output  measurements 
is  currently  under  investigation  at  Hughes  Aircraft.   In  this  method,  dis- 
crete samples  oZ   the  input  and  output  time  histories  are  used  as  coeffici- 
ents in  two  polynomials.   The  output  polynomial  is  periodically  divided  by 
the  input  polynomial  through  the  use  of  a  high  speed  digital  computer.   The 
coefficients  of  the  resultant  polynomial  can  he  shown  to  be  amplitude  sam- 
ples of  the  system  impulse  response.   One  requirement  for  this  method  is 
that  the  amplitudes  of  the  system  output  samples  must  be  corrected  for 
initial  conditions  existing  at  the  beginning  of  the  interval  over  which 
samples  are  taken.   The  correction  process  requires  a  knowledge  of  the 
impulse  response  and  thus  the  impulse  response  must  be  determined  in  a  con- 
vergent succession  of  polynomial  divisions.   At  present,  the  requirements 
for  convergence  to  the  true  impulse  response  are  under  investigation. 

3)  Parameter  Adaptation  Systems 

The  function  of  the  identification  process  in  a  parameter  adaptation 
system  is  the  direct  determination  of  plant  parameters  through  measurements 
of,  and  operations  upon,  selected  signals  and  their  derivatives. 

One  of  the  first  studies  of  parameter  adaptation  was  conducted  by  Kal- 
man,  of  Columbia  University,  who  used  a  digital  computer  to  determine  plant 
parameters  (Ref.  91).   In  Kalman's  system,  coefficients  of  the  plantes  pulse 
transfer  function  were  computed  by  continuously  measuring  and  storing 
sampled  values  of  plant  input  and  output.  After  a  series  of  samples  had 
been  stored,  values  of  the  transfer  function  coefficients  were  assumed  and 
a  series  of  sample  output  points  were  digitally  calculated  using  the  known 
past  input  samples.   The  mean  square  error  between  calculated  and  known  out- 
put points  was  determined  and  the  assumed  values  of  the  coefficients  were 


alt ere  '  ast  squares  filtering  which 

changed  coefficients  in  a  maimer  that  t       the  mean  square  error.   This 
process  was  repeated  at  very  high  speed  until  the  MSE  was  reduced  below  a 
desired  limit.  After  the  plant     Efficients  had  been  determined  with  suf- 
ficient accuracy,  an  optimal  controller  was  synthesized  in  a  second  digital 
computer  routine  and  the  results  were  used  to  adjust  the  system  controller. 
A  different  approach  to  parameter  adaptation  was  taken  by  Whitaker,  of 
MIT,  who  used  a  model  reference  as  shown  in  Figure  26.  (Ref.  171). 
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Figure  26.   MIT  Model  Reference  Parameter  Adaptation  System 
The  input  signal  was  fed  to  both  the  plant  and  to  the  ideal  model,  and  a 
difference  function  between  plant  output  and  model  output  was  used  to  form 
an  integral  squared  error.  The  ISE  thus  formed  was  required  to  be  related 
to  a  single  plant  parameter  which  was  then  varied  in  a  closed  loop  manner 
to  cause  the  tSE  to  attain  its  minimum  value.  The  primary  objective  of  the 
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study  was  the  determination  of  what  functions  of  the  plant-output- to-model- 
output  error  were  indicative  of  what  parameter  variations.   After  the  de- 
sired relationships  were  determined,  the  system  was  mechanised  and  flight 
tested  in  an  F-94A  fighter  aircraft  and  it  was  found  that  adpatation  could 
occur  within  an  interval  corresponding  to  two  or  three  time  constants  of  the 
dominant  system  response  modes.   However,  recent  developments  have  made 
possible  a  significant  reduction  in  convergence  time  and  the  elimination 
of  much  of  the  complexity  of  the  previous  mechanisation  (Ref.  127).   These 
developments  were  derived  from  a  recent  doctoral  research  program  at  MIT 
which  evolved  an  analytical  design  method  for  parameter  adaptation  model 
reference  systems. 

A  method  of  adaptation  employing  an  analog  computer,  which  calculated 
values  of  Instantaneous  forward  loop  gain,  was  studied  by  Corhin  of  Broo!  - 
lyn  Polytechnic  Institute  (Ref.  38),  The  gain  value  calculated  was  com- 
pared with  a  desired  level  and  the  difference  used  as  a  signal  to  the  gain 
controller.  The  differential  equations  of  motion  of  the  plant  were  set  up 
in  the  analog  computer  as  integral  equations  and  were  solved  for  the  gain 

alue,  given  continuous  measurements  of  the  plant  input  and  output  signals, 
as  shown  in  Figure  27. 
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Figure   27. 


Corbin's   Gain  Adaptation  Method  Utilizing  an  Analog  Computer 
in   the   Identification  Process 
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A  method  similar  to  CorMn's  approach,  but  which  employed  a  digital 
computer  to  directly  solve  for  the  plant  parameters,  was  investigated  by 
Azgapetian  of  Servomechanisms  Incorporated  (Ref.  7.).   fri  '"his  study,  the 
equations  of  motion  of  the  second  order  system  under  investigation  were 
converted  to  sample  data  notation,  and  relationships  for  damping  ratio, 
natural  frequency,  am!  forward  loop  gain  were  derived  in  terms  of  past 
sample  points  of  plant  input  and  output.  The  equations  were  then  mechaniz- 
ed in  a  digital  computer  routine  which  periodically  sampled  the  input  and 
output";  signals  and  computed  plant  parameters  at  a  high  rate.   Parameters 
thus  calculated  were  used  in  a  closed  loop  manner  to  vary  a  controller  which 
caused  overall  system  dynamics  to  remain  constant  as  plant  characteristics 
were  varied. 

The  adaptation  method  investigated  in  this  study  is  similar  to  that 
discussed  in  the  preceeding  paragraph  in  that  a  digital  computer  is  util- 
ised to  solve  second  order  equations  of  notion  for  system  damping  ratio 
and  natural  frequency  given  system  output  samples.   The  primary  differences 
are  the  method  of  solution,  and  the  fact  that  the  system  under  investigation 
is  actually  a  sixth  order  system  which  includes  a  second  order  hydraulic 
servo  and  a  fourth  order  airframe  transfer  function. 
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V     0  '  ■   SYSTEM 
As  -mentioned  in  the  preceding  section,  it  is  proposed  that  an  air- 
craft pitch  damping  system  could  be  made  to  produce  invariant  short  period 
damping  ratio  over  all  flight  conditions  by  digital  computation  of  the 
short-period  damping  ratio  and  closed  loop  adjustment  of  the  pitch  damper 
gain  to  cause  the  damping  ratio  to  equal  a  desired  value.  A  block  diagram 
of  the  proposed  system  is  shown  in  Figure  28  below.  The  proposed  identifi- 
cation process  thus  consists  of  solving  the  airframe  equations  of  motion 
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Figure  28.   Proposed  3el£- Adaptive  Fitch  Damper 
for  values  of  damping  ratio  given  samples  of  pitch  angle  and  its  deriva- 
tives, and  the  command  process  consists  of  varying  damper  gain  at  a  rate 
proportional  to  the  error  between  desired  and  computed  /   and  in  the 
direction  indicated  by  the  sign  of  the  error.   Basic  assumptions  and  condi= 
tions  for  this  study  are: 


50 


1,  The  airframe  and  control  system  are  linear.   This  has 
been  shown  to  be  true  for  small  perturbations. 

2.  There  is  no  control  surface  rate  limiting  or  deflection 
limiting. 

3.  Airframe  parameters  are  time  variant  but  are  assumed  con- 
stant over  the  length  of  time  required  for  one  identifica- 
tion cycle. 

4,  System  disturbances  due  to  external  noise  (atmospheric 
turbulence)  will  not  be  considered. 

It  is  also  assumed  that  the  short-period  mode  is  dominant  in  the  iongitudin- 
al  dynamics  and  thus  the  airframe  longitudinal  response  can  be  closely  des- 
cribed by  the  short-period  frequency  and  damping  ratio.   This  assumption 
was  verified  by  conducting  two  longitudinal  flight  simulation  programs  at 
each  of  the  11  conditions  listed  in  Table  II.   The  first  program  was  run 
in  two  degrees  of  freedom,  pitch  and  normal  translation,  and  thus  simulated 
the  short-period  mode  dynamics,  and  the  second  program  simulated  the  short- 
period  and  phugoid  modes  by  inclusion  of  the  longitudinal  translation  degree 
of  freedom.   Damping  ratio  and  period  were  determined  from  pitch  angle  t ' 
histories  obtained  by  pulsing  the  aircraft  about  the  pitch  axis  and  alio1 
it  to  return  to  equilibrium.   Characteristics  obtained  from  the  simulation 
programs  are  listed  in  Table  III.   It  is  seen  that  the  two  degree  and  three 
degree  of  freedom  characteristics  are  in  close  agreement,  and  thus  the  deci- 
sion was  made  to  investigate  a  computational  method  of  determining  the 
!'pseudo"  second  order  damping  ratio  of  a  sixth  order,  three  degree  of  free- 
dom,  airframe-control  system  combination  from  sampled  data. 
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Table  III  -  Short -Period  Mode  Characteristics  Compared  with  Damping  Ratio 
r-,  !  Period  Obtained  frcir.  a  Longitudinal  Three  Degree  of  Free- 
dom F light  S imu ! at ion . 


Fli»h?- 

Gondit ton 

Damping 

Rat  to 

Period, Sec 

"• 

Al t itude 

2  DF 

3DF 

2DF 

3DF 

0,9 

lO.CCO1 

0.444 

'.'■0 

1 .  650 

1.700 

i.2 

10,000 

0.135 

0.177 

1.050 

1 .  005 

1 

10,000 

.  ;55 

0.156 

1.025 

1  020 

r 

25.0001 

245 

0.244 

1.974 

2.( 

1.2 

■t 
25, 

0.100 

0.1028 

1.225 

1.275 

l.l 

25,000 

0.104 

1 .  105 

1.199 

1.175 

1.8 

25,00c1 

0.095 

196 

1 .  200 

1,225 

I  .9 

45, TOO1 

0.140 

' 

2  900 

2.975 

45, C 

0.064 

0.065 

1.725 

1.725 

1.4 

45,  ( 

0.0 

0.067 

1.650 

1.575 

1.8 

45,00c1 

0,  057 

0.058 

1. .  600 

1.599 
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Pre  i  hod 

The  generalized  differential  equation  of  motion  for  a  second  order 
system  is: 

Q  +  LSUnQ+utfGvK^ 

from  which  an  expression  for  danping  ratio  can  be  derived  as: 

r~      K6e  r        _        O          _  CUnB 
2-CUne            ZUJr>e              2.8 
It  is  seen  that,  in  this  form,  the  damping  ratio  is  a  function  of  the  in- 
stantaneous values  of  8  and  its  first  and  second  derivatives,  the  input 
forcing  function,  and  the  natural  frequency.   Azgapetian  has  shown  in  Ref. 
7.   that  an  expression  for,  §       which  is  independent  of  CJn    ,  can  be  de- 
rived as  an  integral  equation.   But  it  was  decided,  for  this  study,  to 
avoid  long  term  integrations  and  their  attendant  drifts,  and  rather  to  make 
use  of  the  high  speed  iterating  capabilities  of  digital  computers  by  de- 
riving a  companion  expression  for  natural  frequency  to  he  used  in  an  itera- 
tive process  for  arriving  at  <f   and  U)n     simultaneously.   The  expression 
for  natural  frequency  being:  i 

which  is  in  terms  of  instantaneous  values  of  8  and  its  first  two  deriva- 
tives, control  input,  and  damping  ratio.   It  is  noted  that  the  term,  k, 
by  which  the  control  deflection  is  multiplied,  corresponds  to  the  airplane 
elevator  effectiveness,  or  steady  state  airframe  gain,  which  is  also  a  func- 
tion of  flight  conditions  and  which  therefore  must  be  determined  in  a 
spearate  computation  routine.   Corbin,  in  Ref.  38.   has  devised  a  method 
of  computing  the  instantaneous  value  of  forward  loop  gain  of  a  second  order 
system  by  utilizing  integral  equations.   Corbin's  system  utilized  an  ana- 
log computer  which  was  supplied  continuous  values  of  plant  input  and  out- 
put, but,  the  method  could  be  converted  to  digital  form  with  sampled  in- 
puts and  used  in  this  study  if  desired.  It  was  decided  not  to  introduce 
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this  additional  complexity  irf:o  the  program  but  father  to  concentrate 
solely  on  the  plant  parameter  computation  method  and  its  application 
to  an  adaptive  system  because  the  computation  method  was,  as  yet,  untried. 
This  decision  necessitated  the  exclusion  of  pilot  control  inputs  to  the 
system  because,  if  pilot  inputs  were  to  be  included,  the  value  of  the  air- 
frame gain  would  have  to  be  known.   It  should  be  made  clear  that  control 
deflections  called  for  by  the  pitch  damper  are  not  included  in  the  preceedinj; 
expressions  for  natural  frequency  and  damping  ratio  because  the  airplane 
pitching  moment  created  by  the  damping  system  is  considered  to  be  added  to 
the  natural  aerodynamic  damping  moment  to  make  up  total  damping.   If  pilot 
inputs  were  to  be  included,  the  control  deflection  created  by  the  pitch 
damper  would  be  subtracted  from  total  deflection  and  the  remaining  control 
deflection  would  be  used  as  the  input  term,  §£   ,  in  the  /  and  d)n      equa- 
tions.  With  pilot  control  input  excluded,  the  equations  are  reduced  to: 


CJ 


«   _  LA-  +     AlL&A2-  -    e. 


f=    - 


e  VI  9  I         ® 

e  OJnB 


The  proposed  iterative  computational  routine  is  as  follows: 

1)  Measure  and  store  sampled  values  of  0,  0,  and  0  at  two 
separate  instants  in  time. 

2)  Assume  a  value  of  damping  ratio  and  compute  ^n      using  the 
first  sample  set  of  0,  0,  &  0  and  assumed  <f   .   The  0Jn 
equation  is  used  in  the  initial  calculation  because  the 
value  of  COr)      is  dependent  upon  variations  in  £  to  a  lesser 
extent  than  is  the  value  of  £  dependent  upon  variations  in 
COn     . 

3)  Use  the  value  of  CJn    computed  in  2)  and  the  second  sample  set 
of  0,  0,  &  0  to  compute   £  . 
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4)  Use  the  ne*  value  of  £  and  the  first  sample  set  to  com- 
pute a  new  value  of  O0n   . 

5)  Repeat  steps  3)  and  4)  in  sequence  until  the  successive  values 
of  damping  ratio  and  natural  frequency  either  converge  or 
diverge. 

Although  it  might  seem  that  only  one  sample  set  of  0  and  its  first  two 
derivatives  is  required  to  define  the  parameters  of  a  second  order  system, 
it  is  seen  that  the  method  described  above  requires  two  sets  for  convergence. 
If  only  one  set  were  used  in  both  equations,  the  OJn    obtained  in  the  first 
calculation,  when  substituted  into  the  second  equation  would  produce  the 
originally  assumed    t,     ,  and  convergence  would  be  impossible.   This  is 
because  the  relationships  for  CJn     and  ^     are  actually  just  two  forms  of 
the  same  equation. 

An  appropriate  question  to  ask  at  this  time  is,  will  the  computations 
converge,  and  if  so,  under  what  conditions  will  they  converge?  The  answer 
is  given  in  the  following  section. 
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VI  APPLICATION  OF  PROPOSED  ITERATIVE  DIGITAL  DAMPING  PATIO  COMPUTER  TO 

A  SECOND  ORDER  SYSTEM 
It  was  felt  that  it  would  be  better  to  test  the  proposed  plant  identi- 
fication method  first  with  a  simple  2nd  order  system  before  adding  the  com- 
plexity of  the  phugoid  and  actuator  modes  of  the  three  degree  of  freedom 
flight  simulation.   When  done  in  this  manner,  the  validity  of  the  basic 
idea  could  either  be  proved  or  disproved  with  certainty  because  the  method 
is  based  on  the  assumption  that  the  airframe  closely  resembles  a  second 
order  system.   Objectives  of  the  program  were  to  determine: 

1)  Will  the  iterative  computing  method  converge  to  correct 
values  of   /  and  cdn     ? 

2)  Conditions  for  convergence. 

3)  Should  both  forms  of  the  0)n    equation  be  used,  and  if  so 
when  is  a  particular  form  applicable? 

4)  Effect  of  system  damping  ratio  on  convergence. 

5)  Effect  of  sample  spacing  on  convergence. 

The  second  order  system  used  in  the  following  programs  was  of  the  form: 

§  «■  (DAMPING    FACTOR ) (I I)  8  +3G 8  =  0 
where  the  natural  frequency  =  6  radians  per  second  was  approximately  the 
median  of  the  airframe  short  period  (*Jn    of  this  study,  and  the  damping 
factor  could  be  varied  to  produce  damping  ratios  over  the  range  from  0.1 
to  i.O  thus  approximating  the  variation  of  the  short-period  mode  damping 
ratio.   Outputs  were  obtained  by  releasing  the  system  from  an  initial  dis- 
placement, 0O  at  time  =  0  and  allowing  it  to  oscillate  freely.   Samples  of 
9,    9,  &   8  were  periodically  taken  and  stored  for  use  in  the  /  calculations. 
The  dynamic  simulation  was  used  in  four  programs  which  were  performed  to 
test  the  damping  ratio  calculator  and  the  influence  of  certain  variables 
on  its  operation.   The  programs  and  results  are  discussed  in  the  following 
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pages  and  sample  listings  of  each  are  included  in  appendix  VI. 

Program  1.   Calculation  of  C0n    Using  Known  Values  of  /  , 
and  Calculation  of  ^   Using  Known  Value  of  tUn  , 
As  a  check  on  the  accuracy  of  the  dynamic  simulation,  the  samples  of 
output  were  used  in  a  series  of  non-iterative  calculations  in  which  U)n 
was  determined  using  the  known  value  of  £  and  ^  was  determined  using  the 
known  value  of  U)n   . 

It:  is  recalled  that  the  relationship  for  natural  frequency, 


«—¥-*  i/fW-4 


e    ^  yv  e  i        e 

has  two  forms  which  will  be  referred  to  for  identification,  as  the  posi- 
tive and  negative  forms.   Both  forms  were  used  in  the  calculations  in  order 
to  determine  under  what  conditions  each  was  applicable.   A  third  purpose 
of  the  program  was  to  obtain  a  basis  for  comparing  the  accuracy  of  f    and 

con      computed  by  the  iterative  method.   Nine  runs  were  made  in  which  damp- 
ing ratio  was  varied  from  0.1  to  0.9  in  0.1  steps,  and  partial  listings 
of  the  runs  made  at  ^  =0.2  and  0.7  are  presented  in  table  IV. 

It  is  seen  from  table  IV  that  the  calculated  values  of  ^   and  OJn 
have  approximately  0.5%  accuracy  except  in  regions  where  the  ratios  -2-  and 

— g—  become  large.   The  errors  in  {      and  U)n      obtained  with  this  series  of 
runs  are  due  to  errors  in  the  generated  output  values;  8,  G,  &  9;  and  thus 
are  indicative  of  the  solution  accuracy  of  the  simulated  second  order  system, 
and  they  also  represent  the  highest  accuracy  which  can  be  expected  from  the 
iterative  method.   It  is  also  observed  that  the  calculated  values  of  £  have 
very  little  variation  whereas  the  cdn    error  increases  in  the  regions  where 
0  is  small,  causing  the  ratios  of  ~-   and  -^—  to  be  large,  and  it  is  in 
these  regions,  before  G  changes  sign,  that  the  negative  form  of  the  uJn 
equation  is  applicable.   Increasing  damping  ratio  had  the  effect  of  widen- 
ing the  regions  in  which  0Jn    errors  appear  and  in  which  the  negative  U)n 
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Table  IV   Calculation  of  Natural  Frequency  Given  Correct 
Damping  Ratio,  and  Calculation  of  Damping  Ratio  Given  Correct 
Natural  Frequency  for  a  Second  Order  System  With  COn-  6.0  Rad/ 
Sec  and  5  =  0.2  and  0.7 

Sample  Spacing  0.0125  Cycle. 


PLANT  IDENTIFICATION  ROUTINE  -  SECOND  ORDER 


CONDI TI0N= 


2  DAMPING  RATI0= 


-2000 


ACCELERATION 

VELOCI  TY 

DISPLACEMENT 

OMEGAN1 

0MEGAN2 

DAM PR AT 

-34.9098  19 

-.434735 

.997322 

6.004203 

-5.829842 

.  1  A0'*94 

*33. 595382 

-.863053 

.989194 

6.004834 

-5.655841 

.  1946  19 

-32. 133230 

-1.274004 

.975818 

6.005478 

-5.483248 

.  195985 

-30.5358  1 1 

-1.665819 

.957424 

6.006144 

-5.310185 

. 19  66  72 

-28.815956 

-2.036888  +itntL   .934259 

6.006840 

-5.  134  753 

.  19  70  89 

-26.986802 

-2.385762 

.906594 

6.007577 

-4.954950 

.19  73  72 

-25.061 701 

-2. 7  11  158 

.874713 

6.008367 

-4.768574 

.197578 

-23.054149 

-3.01 1961 

.838918 

6.009225 

-4.573107 

.197  736 

-20. 977702 

-3.287225 

.79952  1 

6.0101 72 

-4.365574 

.  197863 

-18.845898 

-3.536  1  73 

.756847 

6.01 1234 

-4. 142336 

.  19  79  68 

-16.672 187 

-3.758197 

.71 1229 

6.012447 

-3.8988  10 

.  1980  5  7 

-14.469854 

-3. 952858 

.663006 

6.0  13863 

-3.629053 

.  198134 

-12.251956 

-4.  1  19879 

.612522 

6.015557 

-3.325  122 

.  198203 

-10.031248  - 

-4.259145 

.560125 

6.017644 

-2.976075 

. 19  8264 

-7.820131 

-4. 3  70  700 

.506159 

6.020313 

-2.566302 

.  198521 

-5.630584 

-4.454739   , 
-4.5  11603 

,   .450972 

6.023889 

-2.072655 

.  198  3  73 

-3.474 1 19 

.394904 

6.028996 

-1.459  177 

.198422 

-1.361722 

-4.54  1775 

.338293 

6.036992 

-.666768 

. 19  84  69 

.696184 

-4.545873 

.281469 

6.051495 

.408725 

.198514 

2.689788 

-4.524639 

.224752 

6.086341 

1.966338 

.  198558 

4.609921 

-4.478935 

.  168455 

6.272t>09 

4.3628  38 

. 198602 

6.44  8092 

-4.409  733 

.  1128  7  7 

9. 794 177 

5.832562  • 

.  198645 

.  8.196513 

-4.318107 

.058305 

23.690226 
5.9  72  7  16 

5.934  103 
-4  0.82 1967 

.  198688 

1 1. 596610 

-4.0  72334 

-.046742 

.  198  7  78 

12.836410 

-3.920762 

-.0967  18 

5.9  79  7  35 

-22. 194992 

.  198825 

14. 162731 

-3.751897 

-.  144689 

5.984347 

-16.3  56640 

.  1988  75 

15.371551 

-3.567184 

-,  190449 

5.987638 

-13.479787- 

. 198929 

16.459617 

-3.3681  12 

-.233809 

5.990128  . 

rl 1. 752290 

. 19  89  86 

17.424440 

-3. 156207 

-.274598 

5.992097 

-10.589657 

.199049 

18.264289 

-2.933022 

-.312667 

5.9937D8 

-9.745970 

.199120 

18.978178 

-2.700  125 

-.347884 

5.995064 

-9.099695 

. 199200 

19.565846 

-2.459094 

-.380136 

5.996233  " 

-8.583825 

.  199292 

20.027742 

-2.21 1504 

-.409334 

5.997261 

-8.  158339 

.199402 

20.364998 

-1.958920 

-.435403 

5.998182 

-7.797821 

. 199535 

20.579404 

-1.702891 

-.458292 

5.999021 

-7.485314 

.  199704 

20.673379 

-1.444938 

-.477967 

5.999795 

-7.209030 

.  199725 

20.649939 

-1. 186546 

-.4944  14 

6.000520 

-6.960482 

.200234 

20.512663 

-.929164 

-.507635 

6.001207 

-6.733358 

.200700 

20.265655 

-.674188 

-.517653 

6..001866 

-6.522824 

.201496 

19.913508 

-.  422961 

-.524506 

6.002506 

-6.325066 

.203192 

19.461264 

-.  176  767 

-.528248 

6.003134 

-6. 136986 

.209474 

18.914369 

.063176 

-.528951 

6.003 7  58 j 

-5.955983 

. 168654 

18.278638 

.295721 

-.526699 

6.004383 

-5.779798 

.  192338 

17.560206 

.5  19797 

-.521593 

6.005016 

-5.606394 

. 195 131 

16.765485 

. 734408 

-.513744 

6.005666 

-5.4  3  38  57 

. 196223 

15.901 123 

.938644 

-.503276 

6.006340 

-5.260313 

.  196^5  10 
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Table  IV  -   Continued 


ACCELERATION 

14. 973961 

1 3.990984 

12.959284 

1  1  .886012 

10.778344 

9.643431 

8.488369 

7.320 155 

6.  145654 

4.971565 

3.804387 

2.65039  1 

1.515589 

.405  7  1 1 

-.673819 

-1.717905 

-2.721 799 

-3.681 1 13 

-4.591836 

-5.450344 

-6.253409 

-6.998207 

-7.682317 

-8.303/29 

-8.860837 

-9.352439 

-9. Z77/35 

-10.  136312 

-10.428141 

-10.653565 

-10.813283 

-10.908339 

-10.940106 

-10.910265 

-10.820773 

-10.673937 

-10.472195 

-10.218.301 

-9.915192 

-9.565996 

-9. 174003 

-8.742644 

-8.275467 

-7.7761 16 

-7.248304 

-6.695796 

-6.  122380 

-5.531850 

-4.927986 

-4.314526 

-3.695155 

-3.073483 

-2.453023 

-1.837  183 

-1.229240 

-.632336 


1. 
1. 
1. 
1. 
1. 


131 
312 
481 
63  6 
7/8 
1.905 
2.019 
2.  1  18 
2.202 
2.271 
2.326 
2.366 
2.392 
2.404 
2.403 
2.388 
2.360 
2.  320 
2.268 
2.205 
2.  132 
2.049 
1.957 
1.857 
1.  750 
1.636 
1.516 
1.392 
1.26  3 
1.  132 
.997 
.86  1 
.  725 
.588 
.452 
.318 
.  186 
.056 
-.069 
-.  190 
-.308 
-.420 
-.526 
-.  626 
-.720 
-.807 
-.888 
-.960 
1.026 
1.  084 
1.  134 
1.  176 
1.210 
1.237 
1.256 
1.268 


675 
760 
246 
569 
253 
9  14 
255 
068 
232 
711 
550 
876 
890 
869 
160 
172 
380 
313 
554 
735 
527 
643 
825 
846 
500 
598 
966 
4  33 
636 
007 
771 
944 
327 
699 
619 
418 
199 
831 
053 
857 
025 
043 
441 
795 
726 
900 
033 
888 
274 
048 
1  12 
417 
955 
762 
917 
539 


•time 


SPLACEMENT 

0KEGAN1 

0MEGAN2 

OAMPKAT 

.490325 

6.00/046 

-b. 083842 

.19/  180 

.475034 

6.007797 

-4.902394 

.  19/43/ 

.457558 

6.008604 

-4.7  1369  1 

.  19/62/ 

.438058 

6.009486 

-4.515100 

.197/75 

.416701 

6.010463 

-4.303480 

. 197895 

.393660 

6.01 1564 

-4.074955 

.197995 

.369  1 13 

6.0  12829 

-3.824603 

.198080 

.343239 

6.014315 

-3.545987 

.198155 

.316222 

6.0  16  107 

-3.230430 

.198221 

.288245 

6.018338 

-2.865862 

.198281 

.259490 

6.021224 

-2-434886 

.198336 

.230141 

6.025154 

-1.9  1  1379 

. 198387 

.200378 

6.030897 

-1.254149 

. 198U36 

.  170378 

6.040203 

-.394233 

. 193462 

.  140314 

6.058135 

.792692 

.19852  7 

.1  10354 

6. 107544 

2.548848 

.196571 

.08066  3 

6.565521 

5.  139424 

. 1986  14 

.05  1396 

12. 176150 

5.882218 

. 19865/ 

.022704 

34.02  3  7  96 

5.944415 

. 198/01 

.0052.72 

5.964874  • 

-173.323327 

.198/45 

.032396 

5.975077 

-32.305450 

. 198/91 

.058545 

5.9812  37 

-19.985186 

. 198839 

.083600 

5.985395 

-15.352948 

. 198rt90 

.  10/456 

5.988418 

-12.904145 

. 198945 

.  130016 

5.990736 

-1  1.376240 

. 199004 

.151191 

5.992589 

-10.322459 

. 199069 

. 170707 

5.994 1  18 

-9.5445  12 

. 199142 

.  189095 

5.995415 

-8.940881 

. 199225 

.205/01 

5.996540 

-8.454162 

. 199322 

.2206/8 

5.997534 

-8.04940V 

. 1994  37 

.23399  1 

5.998430 

-7.  /04088 

. 199^/9 

.245615 

5.999248 

-7.402979 

. 199761 

.255536 

6.000007 

-7.  135387 

.200003 

.263748 

6.000720 

-6.39  35  39- 

.200347 

.270257 

6.001398 

-6.67160.3 

.200881 

.275075 

6.002051 

-6.465079 

.201841 

.278226 

6.002687 

-6.270382 

.204  106 

.279742 

6.003513 

-6.084575 

.216423 

^279662 

6.00  39  36 

-5.905169' 

. 184 186 

.278033 

6.004563  ■ 

•  >.-5. 729981 

.193503 

.274909 

6.0052.00 

-5.55  70  17  •-. 

.  195-}30 

.270353 

6.005856 

-5.384384- 

.  196424 

.264432 

6.006538  . 

-5.210202 

. 196932 

.257218 

6.007256 

-5.032525 

.197262 

.248789 

6.008021 

-4.849246 

. 197497 

.239228 

6.008848 

-4.658001 

.197673 

.228621 

6.0097  55 

-4.456031 

.197812 

.217057 

6.010/63 

-4.240008 

.  197925 

.204630 

6.01 19  07 

-4.005/96 

.198020 

. 191432 

6.013228 

-3.748096 

.198102 

. 177561 

6.014792 

-3.4599  17 

.198174 

.1631  12 

6.016693 

-3. 131756 

.  198239 

. 148182 

6.019084. 

-2.750264 

.  198297 

. 132870 

6.0222  18'"' 

-2.295984 

.  198351 

. 117270 

6.026563 

-1.739318 

.  198402 

. 101478 

6.033073 

-1.032846 

.  198449 
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Table  IV  -  Continued 


CONDITION^    7  DAMPI  -!G  KATIO= 


7000 


ACCELERATION 
-32.518630 
-29.005397 
-25.687407 
-22.563123 
-19.6299143 
-16.884t>00 
-14.322/59 
-1 1.940102 
-9. 731415 
-7.691 155 
-5.815440 
-4.092098 
-2.520738 
- 1.092rt01 
.198384 
1.359554 
2.397471 
3.318882 
4. 130484 
4.838395 
5.450624 
5.972048 
6.409392 
6. 768709 
7.055865 
7.276529 
7.436159 
7.539P95 
7.593052 
7.6001 17 
7.565747 
7.474262 
7.389753 

7.096868 

6.9  15523 

6.715228 

6.498949 

6.269440 

6.029253 

5.780741 

5.526069 

5.267215 

5.005985 

4.7440  16 

4.482787 

4.223622 

3.967704 

3.716082 

3.469674 

3.229280 

2.995589 

2.769  182 

2.550546 

2.340076 

2.  138085 

1.944809 


_i-^i  l 


VELOCITY 
-.420471 
-. 804694 
-1. 146321 
-1. 447687 
-1.711196 
-1.939218 
-2.  134075 
-2.298034 
-2.433503 
-2.542022 
-2.626260 
-2.688010 
-2. 729  188 
-2. 75  162  7 
-2.  757078 
-2.  747209 
-2.723603 
-2.687758 
-2. 641089 
-2. 584926 
-2.520520 
-2.449057 
-2.571569 
-2.289128 
-2.202652 
-2. 1 13008 
-2.020993 
■1.927337 
■1. 852705 
1. 757703 
1.642875 
•1.548714 
1.455656 


0 


no  i 


1.  274361 
1. 186762 
1.  101552 
1.  0  18948 
-.9  39  133 
-.862256 
-. 788436 
-.7  17763 
-.650302 
-.586093 
-.525155 
-. 467489 
-.413077 
-.36  1885 
-. 313867 
-.268962 
-.227  100 
-.  188202 
-. 152  180 
-.  1  18940 
-.088382 
-.  060402 
-.034  893 


time 


ISPLACtM 
.997379 
.9896  7  7 
.977440 
.961 187 
.941405 
.918555 
.89306  3 
.865532 
.835752 
.8046  10 
.772284 
.73904  7 
.705169 
.670895 
.63644  9 
.602032 
.567826 
.533993 
.500677 
.468006 
-436089 
.405022 
.374888 
.345754 
.3  176  76 
.29070  1 
.264861 
.240185 
.216682 
. 194566 
. 173238 
. 153292 
. 134516 

_   lUWOh 
.  i  I  <j  kj   /  -r  ' 

. 1004  06 

.085026 

.070727 

.057477 

.0452142 

.035786 

.023672 

.014262 

.005715 

-.002009 

-.008951 

-.015152 

-.020652 

-.025492 

-.029712 

-.033552 

-.056449 

•-059042 

-.04  1  166 

•.042858 

.044  151 

.045078 

.045671 


ENT 


6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

7 

8 

9 

._  in 

i  \; 
11 

13 
15 
18 
23 
29, 
40. 
64. 
153, 

6. 

6. 

6, 

6. 

6, 

6, 

6, 

6. 

6. 

6. 

6. 

6. 

6. 

6. 


0MEGAN1 

.012666 

.0  12680 

.012694 

.012708 

.012722 

.012737 

.012753 

.012769 

.0127  86 

.012804 

.012823 

.012843 

.012865 

.012890 

.012917 

.012947 

.012781 

.015020 

.0  13065 

.013120 

.0  13  186 

.013269 

.013377 

.013524 

.0  13737 

.014078 

.014711 

.016515 

.028057 

.510^41 

.266745 

.155535 

.  138650 

.  ->*. »  -J  <j\i 

.75  7073 

.528653 

.792553 

.807  121 

.049256 

.507034 

.6162-90 

.445604 

.293892 

,0  1240  1 

.012427 

.0.124  50 

.012472 

,0  12492 

,0125  10 

012b27 

012544 

012560. 

0  125  75  , 

012590' 

012604 

012618 

012632 


0MEGAN2 
-5.422600 


-4 
-4 
-5 
-5 
-3 
-2 


8  74  358 


3  7080  3 
904  104 
467957 
05  7  109 
667295 
-2.294832 
-1.9  565  70 


-1 
-1 


6 

6 

6 

6, 

6. 

N  6. 

6. 

6. 

6. 

-414. 

-88. 

-49. 

-54. 

-25. 

-20. 

-17. 

-14. 

-12. 

-1  1. 

-9. 

-8. 

-7. 

-7. 


589  752 

25  1924 

-.920361 

-.594501 

-.270896 

.051857 

.375569 

.702179 

1.0  336  26 

1  .371978 

]  .713973 

2.070580 

2.452076 

2.84  3  154 

5.255442 

3.693355 

4.  162070 

4.667840 

5.217939 

5.813242 

6.006137 

6.0099  10 

6.010925 

6.0  1  139  7 

A  n 1  1  A  70  ■ 

<'  .  v>  i  i  \i  i   ^_ 

.0  1  185-2 
.0  1  1979 
.0  1207  5 
.0  12  150- 
.0122  1  1 
.012261 
.0  12304.. 
.012341 
.012575 
.409029 
.  150044 
.208001 
.015191 
.836909 
.801486 
.  50  269  3 
.7354  17 
.76  1324 
.  187999 
897905 
8  15  159 
888553 
082255 


DAMP.^AT 

.671 594 

.685868 

.690646 

.693036 

.69><'i  78 

.695446 

.696  145 

.6966  76 

.69  7096 

.69 7m 3 7 

.697721 

.697763 

.698  173 

.698  558 

.698522 

.6986  70 

.698805 

.698729 

.699045 

.699  153 

.699254 

.69955  1 

.699444 

.699535 

.6996 19 

.699703 

.699786 

.699367 

.699948 

. 700029 

.700  1  10 

.700192 

.700275 

/nnui 

.   I   VJ  W  ~J   KJ      I 

.700U48 

.700538 

.700635 

.700  73  1 

.700836 

.700947 

.701066 

.70  1 194 

.70  1555 

.70 14  89 

.701662 

.70  1856 

.702079 

.702338 

.702645 

.705017 

.70  34  79 

.704074 

.704872 

.7060  10 

.  70  7  7  73 

.710899 

.718  024 
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Table  IV  -   Continued 


CCE LI-RAT  I  ON   VELOCITY 


1.  76 0'«  12 
1.584997 
1.418607 

1 


261231 

1. 1 12814 

.973255 

.842418 

. 720  133 

.606202 

.500400 

.402483 

.312  188 

.229238 

.153343 

.084206 

.021523 

-.035015 

-.085719 

-. 130900 

-.  170867 

-.205^27 

-.236382 

-.262529 

-.284657 

-.303050 

-.317981 

-.329718 

-.3385  15 

-.344620 

-.348268 

-.349685 

-.349086 

-.346677 

-.342650 

-.337189 

-.330467 

-.322644 

-.3138  74 

-.304298 

-.294047 

-.283244 

-.272002 

-.260425 

-.248609 

-.236640 

-.224599 

-.212556 

-.200577 

-.  188719 

-.  177035 

-.  165568 

-.  154359 

-.  143443 

-.  132848 

-.  122599 

-.  112717 


-.0117 
.  0091 
.02  79 
.0446 
.0594 
.0725 
.0838 
.0936 
.  1018 
.  1087 
.  1  144 
.  1  188 
.  1222 
.  1246 
.  1261 
.  1267 
.  1266 
.  1259 
.  1245 
.  1226 
.  1205 
.  1  175 
.  1  144 
.1110 
.  1073 
.  1034 
.0993 
.0952 
.0909 
.  0866 
.0822 
.0778 
.0  735 
.0692 
.  0649 
.060  7 
.0567 
.0527 
.0488 
.0451 
.0415 
.  0380 
.034  7 
.0315 
.0285 
.0256 
.0228 
.0205 
.0  178 
.0155 
.0134 
.0114 
.0095 
.0078 
.0062 
.004  7 


45 
54 
17 
57 
86 
14 
53 
1  1 
92 
99 
34 
93 
70 
54 
32 
86 
95 
35 
75 
84 
24 
55 
32 
09 
32 
47 
95 
16 
44 
11 
46 
77 
27 
17 
66 
92 
09 
30 
66 
26 
17 
46 
18 
37 
04 
21 
89 
07 
74 
88 
4  7 
48 
87 
60 
64 
94 


time  _' 


? 


LACEMEiMT 

'  0MEGAN1 

0MEGAN2 

OAMPRAT 

045960 

6.0  126  45 

-6.570415 

. 7510  1  1 

045974 

6.012659 

-5-73389  1 

.637835 

045740 

6.012673 

-5.  158  183 

.680^93 

045285 

6.012687 

-4.632083 

.688o02 

044632 

6.012701 

-4.  146765 

.691943 

043805 

6.012715 

-3.695  16  1 

.695794 

042826 

6.0  12  729 

-3.271520   . 

.694977 

04  17  15 

6.  .012  745 

-2.871087 

.695802 

040492 

6.0  12  760 

-2.489869 

.6964  12 

039  174 

6.012777 

-2. 124461 

.696  885 

0377  7  7 

6.012794 

-1.771913 

.697/!65 

036318 

6.012812 

-1.429621 

.697577 

034809 

6.012832 

-1.095246 

.697840 

033265 

6.012853 

-.766647 

.69  80  66 

031697 

6.0  128  77 

-.44 1823 

.698263 

030  1  15 

6.012902 

-.  1  18860 

.698437 

028530 

6.012930 

.204  1  11 

.698594 

026951 

6.012962 

.528959 

.698735 

025384 

6.0129  98 

.85  7599 

.698865 

023838 

6.013040 

1.  192u52 

.698985 

022319 

6.013090 

1.5344  00 

.699096 

020832 

6.013  149 

1.88704  1 

.6992  0  1 

019382 

6.013222 

2.252559 

.69950  1 

017972 

6.0135  16 

2.63  39  04 

.699595 

016608 

6.013440 

3.034477 

.699486 

015290 

6.013614 

3.458259 

.6995  74 

014022 

6.0138  76 

3.9099  75 

.699659 

012806 

6.014320 

4.595298 

.699742 

01 1642 

6.015248 

4.92  10  1  1 

.699-^24 

010532 

6.0184  06 

5.494219 

.699906 

009477 

6. 1588  36 

5.991 120 

.699>86 

0084  76 

6.854  105 

6.008693 

.700067 

007530 

7.6598  54 

6.010503  • 

.700  149 

006638 

8.587596 

6.01  1  184 

.700231 

005799 

9.6  71900 

6.0  1 1543 

.700515 

005013 

10.964613 

6.01  1  7"65 

.700401 

004279 

12.54  1722 

6.011917  « 

.7004  90 

•003595 

14.521378 

*   6.012028 

.70  05  82 

002960 

17.097171 

->  6.0121  13 

.7006  79 

002373 

20.609855 

6.012180 

.700780 

00  1832 

25.7201  14 

6.012236" 

.  700887 

001335 

33.899489 

-   6.0  122  82 

.701002 

000880 

49.225428 

6.012322 

.701  125 

0004  66 

88.733917 

6.012356 

.701259 

000091 

432.9965  7-9 

6.0  12387 

. 7014  06 

000247 

6.0T24  14  - 

-151.100646 

.701568 

000550 

6.012438 

-64.247732 

.701750 

000820 

6.0 12461 

-40.679637 

.701957 

001059 

6.012481 

-29.651905 

. 702196 

001268 

6.012500 

-23.229714 

.702476 

001449 

6.0125  18 

-19.006207 

.7028  1  1 

001604 

6.012535'. 

-16.002548 

.703221 

001736 

6.012552  * 

-13.745592 

. 703739 

001845 

6.012b67 

-1  1.978633 

.704418 

001933 

6.012582 

-10.550309 

.705554 

002002 

6.012596 

-9.365561 

.706735 
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equation  is  applicable.   It  would  be  expected  from  the  results  of  this 
series  of  runs  that  the  iterative  method  of  determining  g       and  Uln    might 
experience  difficulty  in  converging  when  G  is  small  and  that  the  region  of 
non  convergence  will  become  wider  as  system  damping  ratio  is  increased* 

Program  2.   Iterative  Calculation  of  ^  and  (*Jn    Using 
Both  Forms  of  the  CU„   Equation. 

The  purpose  of  this  run  was  to  test  the  iterative  damping  ratio  com- 
puter described  in  section  V.   The  previously  described  second  order  system, 
with  cJn  -  G      &  £=£>2.  ,  was  used  to  supply  output  samples  of  9,  0,  and  © 
which  were  stored,  two  sets  at  a  time,  and  used  as  inputs  to  the     comput- 
er.  Sample  spacing  in  this  run  was  0.05  cycle  and  eacb  $-QJn   calculation 
was  iterated  20  times  using  each  of  the  two  forms  of  the  0Jn    equation.   At 
the  termination  of  each  iterating  process,  the  results  of  each  of  the  20 
iterations  was  printed  out  for  study,  the  oldest  set  of  samples  of  0,  0,  & 
0,  discarded,  a  new  set  generated  and  stored,  and  the  /-^  calculation  cycle 
repeated  using  the  two  new  sample  sets.   A  sample  series  of  the  printouts 
from  the  £-U)n     computer  is  listed  in  Table  V  where  the  results  of  one  con- 
vergent and  one  divergent  solution  are  shown. 

The  primary  result  of  this  program  was  the  demonstration  that  the  itera- 
tive method  would  converge,  and  it  did  so  for  51%  of  the  sample  sets  (^  - 
0.2).   Accuracy  obtained  was  of  the  same  order  as  that  obtained  in  Program  #1, 
approximately  0.05%,  and  in  regions  of  small  values  of  0  the  accuracy  of  the 
iterative  method  was  much  greater.   As  an  example,  COn   from  the  iterative 
method  was  6.00359  as  compared  to  cOn  ~   6.086341  obtained  from  Program  #1 
when  the  identical  sample  set  was  used  with  each  method  in  a  region  of 
small  values  of  0.   It  was  found  that  a  necessary  but  not  sufficient  condi- 
tion for  convergence  of  the  iterative  method  was  that  0  and  0  be  of  opposite 
sign. 


As  is  shown  in  Table  V,  another  interesting  property  of  the  iterative 
method  is  that  when  convergence  occurs,  it  occurs  simultaneously  with  both 
forms  of  the  t*Jn    equation,  and  the  iterative  series  using  the  negative  form 
of  the  o)n    equation  converges  to  the  negative  values  of  $     and  <&#    arrived 
at  using  the  positive  &4j  equation.   Also,  it  was  found  that  the  use  of 
the  negative  ^n    equation  would  not  be  required  in  future  programs  because 
the  iterations  performed  when  using  the  positive  form  converged  in  areas 
where  the  results  of  Program  #1  indicated  that  the  negative  <%  equation  was 
applicable.   It  had  thus  far  been  established  that,  at  £  s  0.2,  the  itera- 
tive computations  would  converge  for  a  length  of  time  corresponding  to  approxi- 
mately one  quarter  of  the  system  period  at  a  rate  of  two  times  per  system 
cycle.   It  was  felt  that  this  would  suffice  for  the  proposed  adaptive 
system,  however,  it  was  considered  desirable  to  determine  the  effect  of 
different  values  of  damping  ratio  on  convergence. 

Program  3.   Iterative  Calculation  of  >  &  ^n   at  9 

Values  of  System  Damping  Ratio. 
Nine  runs  were  made  in  this  program,  with  system  damping  ratio  varying 
from  0.1  to  0.9  in  0.1  steps,  and  a  sampling  interval  =  0.05  period.   Only 
the  positive  form  of  the  cJn    equation  was  used  in  the  iterative  process  and 
the  solution  was  tested  for  divergence  after  each  iteration  to  prevent  un- 
necessary calculations.   Table  VI  is  a  partial  listing  of  the  runs  made  with 
/    =  0.2  and  0.7.   It  is  seen  that  the  region  of  divergence  is  appreciably 
wider  with  £     -   0.7  than  with  0.2.   The  percentage  of  convergent  iterations 
to  total  iterations  is  tabulated  below  as  a  function  of  damping  ratio  for 
each  of  the  nine  conditions; 
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Ratio  yvslltrRllZ%CalCUl$ti:on   of  Neural  Frequency  and  Damping 
&&  ^^/.F^V^F??**-    of  a  Second  OrdePr  6 
Sample  Spacing = 0.05  Cycle 


TIME 
Q.2000G0 

ACCELERATION       VELOCITY            DISPLACEMENT 
-5.630577           -U.45473-..              G.45G971 

MEGAN 

DA.MPR 

7.629993 

0. 303380 

6.818375 

D« 252439 

■-  —  -  v  w    vy     t     V.'  »«/    |     „ . 

6.41318U 

0.22.6C73 

.  6.210091 

3.21257? 

i    **~ 

6.107890 

0. 205695 

i    roy-m 

6.056324 

5.202204 

6.030267 

0.200434 

_  6.017089 

0.199537 

6. OK  U25 

3. 199 : 83 

6.007049 

rJ.  198854 

6.005342 

i.  198737 

6.004477 

9. 198678 

6.004040 

o  1986U8 

6.003819 

).  198633 

60  003707 

0.198626 

6.QC3650 

;. 198622 

6.003621 

&. 198620 

6.0*  36C 7 

D.  198619 
.•.198618 

6.003599 

6.003595 

0.198618 

OH E CAN 

DAMPR 

-2.860373 

0.ft76158 

1  . 

-4.239639 

-0.065534 

-5.07871  ! 

-0.132633 

-5.528898 

-.-.16  55  52       i 

-5.76)842 

-  :o  181965 

-0.190210' 

7-0.T94367 

""  form 

-5. 880 374- 
-5.941390 

-5.972086 

-0.  196468 

-5.987639 

-0.-19753* 

-5.995516 

-0.  19  8.'  67 

-5.999504 

r-'o  198339 

-6.001522 

"''    '   -0.  198477 

-6.002544 

-Q* 198547 

-6.003062 

-0.198582 

-6.003323 

-  )o  1  986  7* ;. 

-6.003456 

-0.  19  86':  9 

-6.003523 

-0.198613 

-6.":'13557 

-0.198616 

-6.003574 

-0.198617 

s 

-6.003583 

-0.198617 

/ 1 
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Table  V   -  Continued 


TIME        ACCELERATION 
0.399999 18.978091+ 

VELOCITY     DISPLACEMENT 
-2-700115    -0.347882 

OMEGAN 
4. 205901 

DAMPR 

0.564624 

3.300829 
2.354981 

0.852037 
1.340582 

1.558809 
0.988132 

2.154067 
3.492864 

0.613964 
0.378345 

5.684419 
9.264260 

0.232391 
0.142551 

15.107418 
24.643762 

(J. 087360 
0.053459 

40.222131 
65.734811 

0.032669 
0.019852 

107.570982 
177.027143 

G. 011444 

0.005859 

307.084564 
599.775604 

0.002441 
-0. 002930 

1439.462204 
-1199.551773 

44936. 139648 
108440.639648 

-2894.778473 
-6985.727661 

261690.726562 

OMEGAN 

-16858.072266 
DAMPR 

-19.193211 

-41.399521 

1.053100 

2.582062 

-97.553423 
-234.412361 

6.248351 
15.085813 

-565.269005 
-1363.942886 

36.408317 
87.862394 

-3291.417694 
-7942.874573 

212.031506   , 
511.678204  r 

-19167.864990 
-46256.204590 

1234.790436 
2979.816833 

-111626.239258 

-269378.281250 

7190.943481 
17353.303223 

-650068.  140625. 
-1568755. 156250 

41877.279297 
1C 1058.941406 

-3785745.781253 
-9135824.375000 

243877.103516 
588528.257812 

-2204672C.5fn0QO   . 
■!4!6lifttf66l!flft9 

-128391576.0r.n00O 
-77459164.000^00 

1420246.000000 
34373<S0«906P54» 
8270963.625000 

19959625.000000 
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Table  VT.  -  -nt  Iterations  as  a  Function  of  System 

Damping  Ratio 

Damping  Ratio  0.1   0.2   0.3   0.4   0.5   0.6   0.7   0.8   0.9 
%  Convergence  52%   517.   51%   49%   48%   41%   34%   27%   19% 
The  percentage  of  convergent  iterations  drops  off  sharply  for  values  of 
exceeding  0.5  but  it  is  not  desired  to  increase  the  airframe  damping  ratio 
above  0.7  and  thus  the  minimum  obtainable  ratio  of  convergent  to  total  itera- 
tions should  never  be  less  than  approximately  1/3. 

Program  4.   Influence  of  Sample  Spacing  on  the  Percentage 

of  Convergent  Iterations. 

The  final  program  was  run  to  determine  what  influence^  if  any,  the  width 

of  sample  spacing  has  on  the  ability  of  the  iterative  calculations  to  converge. 

All  runs  were  made  with  £    =  0.2,  and  six  values  of  sample  spacing  ranging 

from  0.01  cycle  to  0.5  cycle  were  investigated.   Results  of  the  program  are 

tabulated  below. 

Table  VIII  -  Percent  of  Convergent  Iterative  Calculations  as  a  Function  of 
Sample  Spacing  /  =  0.2 

Spacing        0.01  cycle     0.02     0.05     0.1     0.2     0.5 
%  Convergence  56%  53%      51%      ^  1%         48%     46% 

It  is  seen  that  narrow  sample  spacing  favors  convergence s    however  it  can  be 
varied  over  wide  limits  with  little  effect  on  convergence.   It  was  arbitrar- 
ily decided  to  set  sample  spacing  at  0.025  seconds  for  the  three  degree  of 
freedom  airframe  simulation;  the  primary  consideration  being  conservation  of 
computer  time.   This  sample  spacing  corresponded  to  approximately  50  samples 
per  short  period  cycle,  depending  on  flight  condition. 
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Table  VII  Iterative  Calculation  of  Natural  Frequency  and  Damping 
Ratio  of  a  Second  Order  System  With  «4>=6.0  Rad/Sec  and  f-   0.2  &  0.7 
Sample  Spacing  =  0.05  Cycle 

Condition  2    Damping  Ratio  »  0.2 

TINE  ACCELERATION   VELOCITY     DISPLACEMENT 

C. 200000  -5.630577     -4.454730      0.450971 

ONEGA-N  DAMPRAT   „__ 

7.629993  0.303380 

6.818375 0.252439 

6.413184  0.226073 

6.210091  0.212571 

6.107890  0. 205695 

6.056324  0.2C2204 

6.030267  0.200434 

6.017089 0.199537 

6.010423  0.199083 

6.007049  0.198854 

6.005342  0.198737 

6.004477  0.198678 

6.004040  0.198648 

6.003819  0.198633    

6.003707  0.198626 

6.003650  0.198622 

6.003621  0.198620 

6.003607  0.198619 

6.003599  0.198618 

6.003595  C. 198618 


TINE        ACCELERATION   VELOCITY     DISPLACEMENT 

0.250000  2.689775     -4.524627      0.224752  .  • 

0MEGA-N  DAMPRAT 

9.251958  0.261913 i .  '. 

6.966586  0.215692  ^ 

6.256753  0.202902  _z^l 

6.066626  0.199669 

6.019025  0.198874     L_ 

6.007354  0.198680 

6.004508  0.198633 „  .  - 

6.003815  0.198622 

6.003646  0.198619         

6.003605  0.198618 

6.003595  0.198618 

6.003592  0.198618   ~ 

6.003592  0.198618 vf 

6.003592  0.198618 

6.003592  0.198618 
6.0C3592   '   0.198618 

6.003592   j  0.198618 

6.003592   '  0.198618 
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Table  VII-  Continued 


Damping  Ratio  -  0.2 


6.003592 
6.003592 


0.198618 
0.198618 


TIME 
i.  0.299999 

OMEGA-N 

15.324385 

-7.279872 

_-<4. 234596 

-1.191589 


ACCELERATION 

9.848088 

DAMPRAT 

0.085544 

-0.165185 

_-0.  279041 

-0.983381 


VELOCITY 
-H.205210 


DISPLACEMENT 
0.005013 


SOLUTION  DIVERGING 


a  — 

TIME 
J). 349999 
OMEGA-N 


ACCELERATION 
15.371489 


VELOCITY 
-3.567171 


SOLUTION 


DAMPRAT 
DIVERGING 


DISPLACEMENT 
-0.190448 


TIME 
_0.399999_ 
OMEGA-N 
4.205901 
3.300829 


ACCELERATION 
_18._97809Jt 

DAMPRAT 
_„0.  564624 

0.852037 


VELOCITY 
-2.7QQ115 


^SOLUTION  DIVERGING 


DISPLACEMENT 
-0.347882 


TIME 
0.449999 
OMEGA-N 
4.907358 
4.181908 


ACCELERATION 
20.579304 

DAMPRAT 
.  0.570963^ 

0.882182 


VELOCITY 
j=0. 702884 


-SOLUTION  DIVERGING 


DISPLACEMENT 
-0.458290 


TIME 
0.499998 
OMEGA-N 
.5.377643 


ACCELERATION 

20.265548 

DAMPRAT 
__Q.i3iL322 


VELOCITY 
-0. 6741 85 


SOLUTION  DIVERGING 


DISPLACEMENT 
-0.517650 


TIME 


___  ACCELERATION VELOCITY_ 

0.549998     18.278530      0.295719 
_QHE  G  A-N __DARP_RAT 


5.757618    -0.240352 
SOLUTION  DIVERGING 


DISPLACEMENT. 
-0.526696 
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Table  VII-  Continued     Damping  Ratio  ■  0.2 


TIPE        ACCELERATION   VELOCITY     DISPLACEMENT 

0.599997     14.973863      1.131666    -0.490321 

QNEGA-N       DAMPRAT 
6.119884     0.244751 
6.030041      0.209181 
6.009638      0.201036„_ 
6.004975      0.199172 

6.  003908 0.  1  98745   

6.003664      0.198647 

6.  003608____G.  198625 

6.003595      0.198619 

6.003593 0.198618 

6.003592  0.198618 
6.003592  0.198618 
6.003592  0.198618 
6.003592      0.198618 

6.003592 0.198618 

6.003592  0.198618 
6.003592  0.198618 
6.003592  0.198618 
6.003592  0.198618 
6.003592  0.198618 
6.003592 0.198618 


TIME  ACCELERATION   VELOCITY     DISPLACEMENT 

0.649996  10.778271      1.778240    -0.416698 

OMEGA-N  DAMPRAT 

6.52  5987 0.300233 

6.262415  0.249806 

6.132749  0.224382 

6.068288  0.211580 

6.036061  0.205138 : V 

6.019904  0.201897 

6.011790  0.200267 _ 

6.007714  0.199447 

6.005664  0.199035 ' 

6.004634  C. 198828 

6.004116  0.198723 

6.003855  0.198671            .-• 

6.00  3724 0.198645 

6.003658  0.198631 

6.0C3625 0.198625 . 

6.003608  0.198621 

6.003600  0.19862Q__ * 

6.003596  0.198619 

6.003594  0.198618 

6.003593  0.198618 
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Table  VII-  Continued 

Condition  7   Damping  Ratio  =  0.7 


^ACCELERATION 
'  1.359547 

TINE 
0.200000 

VELOCITY 
-2.747204 

DISPLACEMENT 
0.602031 

ONEGA-N 
4.221364 

DAMPRAT 
0.521158 

4.917109 
5.325055 

0.589098 

0.629943 

5.575466 
5.732646 

0.655293 
"'.671299 

5.832533 
5.896477 

C. 68  150b 
0.688051 

5.937594 
5.96U107 

0.692266 
0. 69U986 

5.981233 
„5. 992308 

0.696743 
C. 697880 

5.999476 
6.004117 

0.698616 
0.699093 

6.0C7122 
6.009069 

0.699402 
0.6996C2 

6.C1C331 
6.0111U8 

0.699731 
6.699815 

6.011677 
6.012021 

0. 699870 
C. 699905 

6.012243 

C. 699928 

TINE 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.250C00 

ONEGA-N 

14.  838883 
DANPRAT 

-^.584920 

0.468005 

2.8613U7 
U. 887297 

{..5  86  139 
C. 6  33940 

5.364674 
5.622871 

0.660114 

_T_.675476 

5.773553 
5.864458 

0.684772 
0.690488 

t 

5.920257 
J5. 954843 

0.694035 
C. 696248 

5.976406 
5.989896 

0.697633 
0.698501 ' 

5.998354 
6.003663 

0.699046' 
0.699389 

6.007000 
6.0C9C97 

0.699605 
0.699740 

6.010416 
6.01124  5 

0.69982  5 
U. 69987V 

6.011767 
6.012096 

0.699913 

0.699934 
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Table  VII-  Continued     Damping  Ratio  =  0.7 
6.012302      ;,, 699947 


6.012432      •:«699956 


TINE       -ACCELERATION   VELOCITY     DISPLACEMENT 


0.299999      6.763687     -?. 289121      3.345753 

ONEGA-N       CAMPRAT 
-3.24  9:22     -0.700412 
-1.71765^     -0.990456 

SCLLTICN  CIVERGING 


TIME        ACCELERATION   VELOCITY    DISPLACEMENT 
0.349999      7.539968     -1.927330      5.240162 

OMEGA-N       CAMPRAT 


-9.915033     -0.815083 
„  -2.3*6982     -C. 991636 
SCLLTICN  CIVERGING 


TIME        ACCELERATION   VELOCITY     DISPLACEMENT 
0.399999      7.494233     -1.548708      5.153291 

_  0  f.  E.GA- N CAMPER  A I 


SCLLTION    DIVERGING 


TIME  _A_CCE1.JE.ELAJ  LOIS       VELOCITY     DISPLACEMENT 

0.449999  6.915493  -1.186757  0.085026 

0  MEG  A_-N DAMRPxAX. 

SOLLTICN  DIVERGING 


TIME  .  ACCELERATION   VELOCITY     DISPLACEMENT 

0.499998      6.029223     -0.862252      0.033986  ; 

OMEGA-N     -   CAMRRiLt : '. __ 

-44.580629     -  .957'       .              -  % 
SCLLTICN  CIVERGING : 


TIME        ACCELERATION   VELOCITY     DISPLACEMENT 
0.549993      5.0(5958     -*. 586095  '  '  -!).0Q2„C9 


CMEGA-N       CAMPRAT 
-64.266227      3-0437Q1 


SCLLTICN  CIVERGING 


TIME   ACCE LERATIQN   VELOCITY     DISPLACEMENT 

0.599997      3.967681     -0.361883     -0.G25492 

OMEGA-N  ;  CAMPRAT 

10.228293      0.175711 

20.289338  :    -0.444423 

SCLLTICN  CIVERGING 
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Table  VII-  Continued      Damping  Ratio  ■»  0.7 


TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.649996 

OVEGA-N 

•"  2.995569 
CAMPRAT 

-°.1S82CC 

-0.039041 

8.C2886C 
_   10.427491 

0.158457 
-  }.  318J48 

SOLUTION  DIVERGING 


TIKE-        ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.699996      2.138C70 

-n. 060402 

-0.045078 

OMEGA-N       DAMPRAT 
7.040965     -0.113642 

SOLUTION  DIVERGING 

TIME        ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.749995      1.418596 
ONEGArJ^      camprat 

*. 027917 

-0.045740 

SOLUTION  DIVERGING 

TINE        ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.799994      r.842411 
ONEGA-N       CAMPRAT 

^.083853 

-0.042826 

5.818552      C. 622535 
5.961974 0.679927 


5.999497      0.694778 

6.CC9243  0.698625 

6.011769      C. 699621 

6.012424  0.69988C 

6.012594  <^. 699946 
6.^12637  0.699964 
6.012649       <,  599966 

6.012652  .0.699969 

6.012653  0.69997C 

6.012653 Q. 699970 

6.012653      C. 699970 

_6  .  Q  1  2  6  5  3 _0  .  6  99  9  7  Q 

6.012653      0.699970- 
6. 012653      0.69997,- 
6.012653      0.699970 
6. Q  12653  Q.699'97) 


6.012653      0.699970 
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VII  -  APPLICATION  OF  THE  DAMPING  RATIO  COMPUTER  TO  A  SELF  ADAPTIVE  PITCH 

DAMPING  SYSTEM 

Investigation  of  the  Ability  of  the  ^-^n   Computer  to  Determine 
Short-Period  Mode  Characteristics 

It  has  thus  far  been  established  that  the  damping  ratio  computer  opei- 
ates  satisfactorily  when  supplied  input  samples  generated  by  a  second  order 
system,  and  it  has  also  been  established  that  the  short  period  mode,  a 
second  order  oscillation,  is  the  dominant  mode  of  the  airplane  longitudinal 
dynamics.   The  next  step  in  the  program  is  then  to  investigate  the  opera- 
tion of  the  computer  when  supplied  samples  generated  by  the  three  degree  of 
freedom,  fourth  order  airframe  simulation.   The  simulation  program  is  dis- 
cussed in  detail  in  Appendix  VI. 

Flight  condition  number  seven,  M  =  1,8  @  25,000  ft.,  having  short-period 
mode  damping  ratio  =  0.095  and  period  =  1.2  seconds,  was  chosen  for  the  test, 
At  the  start  of  the  computer  run,  the  airframe  was  released  from  a  non-trim 
pitch  condition  and  allowed  to  oscillate.   Periodic  samples  of  pitch  angle, 
8,  and  its  first  two  derivatives  were  taken  and  supplied  to  the  £-<*4?  com- 
puter as  described  in  the  preceeding  section.   At  the  beginning  of  each  set 
of  £-&n   calculation  an  initial  value  of  f  =  0.4  was  introduced  into  the 
equation  and  the  iterations  were  carried  out  15  times.   The  iterations  con- 
verged for  49  percent  of  the  sample  sets,  which  is  in  good  agreement  with 
the  previous  second  order  system  investigation,  however  the  computed  values 
of  damping  ratio  ranged  from  0.1448  to  0.024565  during  a  2.5  second  (flight 
time)  run.   Such  a  wide  variation  of  computed  £  was  inconsistent  with  the 
results  of  the  preceeding  study  in  which  it  was  demonstrated  that  the  longi- 
tudinal aircraft  dynamics  very  closely  approximate  the  dynamics  of  a  second 
order  system.   A  review  of  foregoing  investigations  was  conducted  to  deter- 
mine if  there  were  any  significant  differences  between  the  second  order 
study  and  the  present  application  of  the  ^-^  computer,  and  the  answer  was 
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found  in  examination  'ar.s  short-period  pitch-angle-to-elevator 

deflection  transfer  fund 

It  is  seen  that  the  1/S  term  causes  pitch  angle  output  to  exhibit  a  theo- 
retically infinite  steady  state  error  or,  if  the  airframe  is  perturbed  from 
its  trim  position  at  a  given  flight  condition^  the  ensuing  short  period  oscil- 
lations will  not  decay  to  the  original  trim  angle.   Furthers  it  can  be  showr. 
that  the  airframe  will  return  to  trim  condition  as  a  function  of  the  longi- 
tudinal phugoid  mode.   The  primary  difference  between  the  ideal  second  order 
system  and  the  short  period  dynamics  is  then  the  fact  that  the  former  system 
oscillates  about  G  equal  zero  whereas  the  airframe  oscillates  about  an  aver- 
age pitch  angle  which  is  constantly  varying  as  a  function  of  the  phugoid 
dynamics,  and  thus  the  pitch  angle  samples  introduced  into  the  computer  were 
erroneous.   The  solution  to  this  problem  was  seen  in  the  pitch-rate-per-stabi- 

lizer-def lection  transfer  function, 

-JL-         *e(6-t-T&) 

which  exhibits  zero  steady  state  error,  and  which  also  shows  that  free  air- 
frame short  period  oscillations  oscillate  about,  and  decay  to,  zero  pitch 
rate,   It  was  realized  that  pitch  rate  is  the  true  second  order  variable  of 
the  airframe  longitudinal  dynamics  and  in  the  study  to  be  outlined  next, 
pitch  rate  and  its  first  two  derivatives  9  and  0  were  utilized  as  inputs  to 
the  %-tOn   computer. 

A  second  test  of  the  $-uin  computer  using  the  three  degree  of  freedom 
flight  simulation  was  performed,  as  stated  above,  and  a  partial  listing  of: 
the  results  obtained  when  sampled  values  of  8,  8  &  8,  were  used  as  computer 
inputs  is  presented  in  Table  IX.   As  before,  the  calculations  converged  for 
49?{,  of  the  sample  sets  but,  as  can  be  seen,  the  calculated  values  of  damp- 
ing ratio  compare  favorably  with  the  known  £  =  0.095.   It  can  also  be  seen 
from  examination  of  Table  IX  that,  in  most  cases,  the  calculations  have  not 


quite  settled  out  to  steady  values  o   £  and  CUn    ai     15  iterations,  and 
as  a  result,  it  was  decided  to  repeat  the  iterative  calculation  cycle  20 
times  in  all  future  programs. 

Application  of  the  /-^  computer  to  the  Adaptation  Loop 

It  was  felt,  at  this  times  that  the  £-^/>  computer  was  generating  suf- 
ficiently accurate  values  of  damping  ratio  to  allow  its  use  in  the  adapta- 
tion loop,  however  one  obvious  discrepancy  remained  to  be  corrected  before 
the  loop  could  be  closed.  That  iSj,  what  should  be  done  regarding  the 
calculations  which  do  not  converge  and  thereby  generate  highly  inaccurate 
values  of  £  ?   Before  discussing  the  technique  which  was  used  to  detect  and 
account  for  divergent  calculations,  possibly  a  better  understanding  may  be 
obtained  from  a  short  review  of  the  proposed  mechanization  of  the  adapta- 
tion loop. 

It  will  be  recalled  that  it  was  proposed  to  compare  the  output  of  the 
damping  ratio  computer  with  the  desired  value  of  short  period  damping  ratio 
•  the  error  between  £  -  desired  and  £  -    computer  being  used  as  the  input 
signal  to  a  rate  servo  which  would  vary  pitch  damper  gain  in  the  proper 
direction  to  null  out  the  £-  error,  and  at  a  rate  proportional  to  the 
magnitude  of  the  error.   It  is  apparent  that  the  erroneous  values  of  £ 
resulting  from  divergent  calculations  could  not  be  allowed  to  reach  the 
computer  output,  and  it  was  found  that  this  could  be  accomplished  by 
simply  including  a  routine  which  examined  the  resultant  values  of  ^  and  £ 
at  the  termination  of  the  series  of  iterative  calculations  for  each  sample 
set.   After  examination  of  a  large  number  of  past  output  records  from  both, 
the  second  order,  and  the  airframe  studies,  it  was  determined  that,  in 
every  case,  when  the  calculations  diverged  the  resultant  calculated  values 
of  either  £  or  COn   or  both  were  completely  out  of  the  range  of  variation 
of  the  real  values  of  £   and  U)n   .   For  examples  of  divergent  calculations, 
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Table  IX  Iterative  Calculation  of  Airframe  Short  Period  Damping 
??  *?  a?d^atUr?;1  Frec^uency  Using  0,  6,  &  0,  Generated  by  a  Long- 
Inputs         Degree  of  Freedom  Flight  Simulation,  as  Computer   S 
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Table  IX  -  Continued 
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see  Tables  V.  or  VII.   Also,  it  was  knovn  with  certainty,  that  for  the  air- 
craft considered  in  this  study,  damping  ratio  vras  restricted  well  within 
the  range  from  C.  to  1,0  and  natural  frequency  within  the  range  from  0. 
to  8.0  radians /second.   Thus  a  divergence  detection  routine  was  programmed 
to  use  a  simple  logical  process  for  determining  whether  or  not  the  calculated 
values  of  both  £  and  OJn   fell  within  the  allowable  ranges.   If  both  values 
were  acceptable  to  the  detection  routine,  the  calculated  value  of  t,    was 
set  into  the  damping  ratio  comparison  process,  and  if  either  £  or  Q)n   or 
both  were  unacceptable  to  the  detection  routine,  the  last  acceptable  com- 
puted value  of  %     was  reset  into  the  comparator.   The  values  of  (x)n   and  £ 
resulting  from  a  divergent  series  were  always  orders  of  magnitude  different 
from  the  allowable  ranges  of  uJn    and   £  ,  and  thus  the  allowable  ranges 
could  be  altered  considerably  to  suit  the  characteristics  of  other  systems 
with  no  deterioration  in  the  effectiveness  of  the  divergence  detection 
routine.   As  an  example,  if  it  were  expected  that  the  airframe  might  have 
short-period  poles  in  the  right  half  plane  over  some  flight  conditions,  the 
allowable  £  range  could  be  increased  from  £  =  -1.0  to  +1.0. 

A  variable  which  has  not  yet  been  discussed,  but  which  was  quite  im- 
portant in  the  evaluation  of  the  adaptation  loop,  is  the  gain  of  the  rate 
servo  (hereinafter  called  the  adaptation  servo)  used  to  vary  pitch  damper 
gain.   When  £-  error  signals  are  supplied  to  the  adaptation  servo  which 
in  turn  causes  variations  in  airframe  damping  ratio,  the  entire  adaptation 
loop  appears  to  resemble  a  second  order  ser*vomechanism  and  should  be  amenable 
to  conventional  systems  analysis  techniques.   However  analysis  was  compli- 
cated by  the  intermittent  action  of  the  damping  ratio  computer  and  thus  an 
"optimum"  value  of  adaptation-sei vo-gain,  K  ,  was  determined  experimentally 
by  performing  a  series  of  three  degree  of  freedom  computer  runs  with  vary- 
ing values  of  K  .   Flight  condition  7,  M  =  1.8  @  25,000  ft.,  was  again  used 
s 
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Figure  29.   Time  History  of  Airplane  &  Control  Surface  Motion 
and  Pitch  Damper  Gain  During  Adaptation  Process  With  No  Adapt- 
ation Servo  Lock 

Initial  K(j  »  0,   Adaptation  servo  gain  =  1. 

Flight  Condition  7  -  M-1.8  @  25,000  ft. 


3      4 
Time,  Seconds 
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for  all  runs  and  adaptation-servo  gain  was  varied  from  K  ~  0.05  to  2.0 

s 

in  six  steps. 

Figure  29  is  a  typical  time  history  of  the  run  made  with  K  -  1.0. 
Airplane  pitch  angle,  pitch  rate,  control  surface  position,  and  pitch  damp- 
er gain  are  shown  as  functions  of  time,  and  again  the  run  was  started  by  re- 
leasing the  airframe  from  a  non-trim  pitch  position  and  allowing  it  to  oscil- 
late.  The  damping  ratio  of  the  basic  airplane  at  this  condition  was  0.096 
and  the  steady  state  damper  gain  for  correct  adaptation^  from  the  root  locus 
plot,  was  0.344.   Initial  damper  gain  was  set  at  zero  at  time  ~   zero  and  it 

is  seen  that  K,  =  0.344  was  achieved  in  0.45  seconds,  and  that  the  airframe 
a 

motion  was  heavily  damped  by  the  control  action.   However,  it  is  also  seen 

that  damper  gain  did  not  level  out  at  the  desired  value  of  K,  =  0.344  but  in- 

d 

stead  continued  increasing  in  a  somewhat  random  manner  until  at  time  s  10 
seconds  (not  shown)  it  reached  the  required  value  for  instability  of  the 
actuator  mode  and  a  violent  divergent  oscillation  ensued.   In  this  condition, 
the  actuator  motions  comprised  the  dominant  longitudinal  mode,  and  it  was 
interesting  to  note  that  the  $~Wn    computer  was  still  performing  its  duty  in 
that  it  was  nox*  computing  natural  frequency  and  damping  ratio  of  the  actua- 
tor mode.   The  crossover  frequency  for  this  mode  was  predicted  to  be  34 
radians  per  second  from  the  root  locus  plot,  x^hich  was  verified  by  the  fact 
that  computed  damping  ratio  changed  from  positive  to  negative  sign  when 
computed  con    was  approximately  34.5  rad/sec.   But,  the  operation  of  the 
adaptation  loop  was  obviously  unacceptable  with  K,  =  1.0,   An  improvement 
in  performance  was  detected  when  lower  values  of  adaptation-servo-gain  were 
investigated,  and  with  K4  =0.1  the  damper  gain  did  settle  out  reasonably 
well,  but  approximately  eight  seconds  were  required  to  complete  the  adapta- 
tion process  when  using  the  K,.   Osburn  has  shown  in  Ref.  12?  that  the  MIT 
self  adaptive  system  is  capable  of  completing  the  adaptation  process  on  a 
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Liar  airfr  two  cycles  of  plant  oscillation,  which  would 

correspond  to  adaptation  time  -   2,4  seconds  for  flight  condition  7, 
presently  umter  liscussion.   Thus,  a  review  was  made  of  the  adaptation 
method  to  determine  if  the  process  could  be  speeded  up  without  attend- 
ant drift  in  pitch  damper  gain. 

The  inability  of  the  adaptation  loop  to  set  proper  damper  gain  when 
operating  at  high  adaptation-servo-gains  was  found  to  be  caused  by  four 
interacting  factors: 

1)  The  reader  will  recall  that  when  the  f-^n   process  failed  to 
converge,  the  computer  set  the  most  recently  computed  con- 
vergent value  of  £  into  the  damping  ratio  comparator,  and 
the  adaptation  servo  continued  to  change  damper  gain  at  a 
rate  proportional  to  the  indicated  £  -  error. 

2)  It  was  found  that  during  the  adaptation  process,  the  const- 
antly changing  damper  gain  caused  errors  in  computed  °^n 
and  £  due  to  the  fact  that  system  parameters  were  altered 
appreciably  during  one  sampling  interval.   And,  the  comput- 
ed values  of  4    were  caused  to  vary  in  a  random  manner  and 
were  generally  0.05  to  0.1  below  the  true  value  of  damping 
ratio. 

3}   As  pitch  damper  gain  increased,  the  effect  of  the  time  rata 
of  change  of  system  parameters  on  the  ^-(x]n     calculations 
was  magnified  due  to  the  increased  system  loop  gain, 
when  system  motions  damped  to  zero,  the  §-(*)„    computer  was 
not  supplied  with  sufficient  information  to  perform  calcu- 
lations. 
It  can  be  seen  how  these  factors  combined  to  cause  erratic  operation  of  the 
adaptation  mechanism.   As  an  example,  if  a  value  of  £    =  0.6  were  computed 
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and  set  into  the  £  comparator,  and  no  new  convergent  values  of  <£  were 
generated  for  tvo  seconds,,  a  steady  £  -  error  signal  =  0.1  would  be  fed 
to  She  adaptation  servo  and  it  would  increase  damper  gain  at  a  constant 
rate  during  the  two  second  period,,   This  is  illustrated  in  Figure  29.   where 

K,  was  increased  linearly  from  time  =  3.5  seconds  to  5,5  seconds 0   At  5.5 

d 

seconds  a  convergent  value  of  £  ,  exceeding  C„7,  was  obtained,  set  into  the 

£  -  comparator,  and  damper  gain  began  a  downward  drift.   However,  the  drift 
of  Kd  was  generally  upward  because  the  erroneous  damping  ratio  solutions 
were  usually  smaller  than  0.7. 

The  problem  of  drifting  damper  gain  was  solved  by  an  alteration  of  the 
%-<*)n    computer  output  routine  in  which  output  logic  was  re-programmed  to 
set  £  =   0.7  into  the  £    comparator  if  no  convergent  solution  were  obtained. 
This  action  effectively  "locked"  the  adaptation  servo  at  its  present  set- 
ting, and  damper  gain  was  not  allowed  to  vary  until  generation  of  a  new 
convergent  solution.   The  above  process  will  be  referred  to  as  the  "adapta- 
tion-servo lock"  in  the  remainder  of  this  report.   An  important  benefit 
accruing  from  the  use  of  the  servo  lock  vac  the  fact  that  the  plant  was 
allowed  to  settle  to  a  linear  system  during  each  period  of  nonconvergent 
calculations  and  thus  the  accuracy  of  subsequent  convergent  calculations 
was  improved.   It  is  recalled  that,  when  the  plant  is  in  motions  the  calcu- 
lations converge  during  one  quarter  period  in  each  half  system  cycle,  and 
from  this,  it  would  be  expected  that  damper  gain  would  be  changed  in  quart- 
er period  steps  during  the  adaptation  process.   This  is  borne  out  in  Figure 
30  which  is  a  time  history  of  the  adaptation  process  with  the  servo  lock  in 
operation. 

The  run  shown  in  Figure  30  is  identical  in  all  other  respects,  except 
the  use  of  the  servo  lock,  to  the  run  illustrated  in  Figure  29  that  is, 
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tion  7  ;  Lon- servo  gain  =1,0.   It  is  seen  that 

pitch  damper  gain  increased  in  steps,  as  predicted,  and  it  settled  out 
close  to  the  predicte    Lu<     EC ,  =  0.344  in  2.3  seconds.   After  three 
seconds  running  time,  the  airplane  motion  damped  cut  and  no  values  of  £ 
were  available  until  t  =  5  sec.   During  this  interval  damper  gain  remained 
constant  at  If,  =  0.3.   At  t  -  5  sec,  the  airframe  was  perturbed  by  an  im- 
pulse and  the  computer  began  generating   values  causing  damper  gain  to  be 
moved  into  the  If,  =  0.34  to  0.35  region  again.   The  time  history  of  motions 
resulting  from  the  second  impulse  also  serves  to  illustrate  the  behavior  of 
the  adapted  airframe  which  displays  a  dan  ping  ratio  =  0.7. 

A  series  of  runs  was  male  with  the  servo  lock  in  operation  using  eight 
different  values  of  adaptation  servo  gain  ranging  from  Ks~   0.1  to  5.0.   It 

was  found  from  this  series  that  values  of  Tf  exceeding  1.0  caused  an  ex- 

s 

cessive  overshoot  in  damper  gain  during  the  adaptation  process  and  that 
values  on  K  less  than  1.0  resulted  in  a  smooth  adaptation  transition,  but 
time  required  to  fully  adapt  was  approximately  inversely  proportional  to 
se'rvo  gain.   On  the  basis  of  these  results,  it  was  decided  that  K  s  1.0 
was  optimum,  and  a  second  series  of  runs  was  conducted  to  investigate  the 
operation  of  the  adaptation  loop,  with  this  gain  setting,  at  each  of  the 
11  flight  conditions.   Table  X  is  a  sample  listing  of  the  computer  output  for 
the  three  degrees  of  freedom  run  performed  for  flight  condition  7,  from  which 
the  reader  can  see  the  stepwise  manner  in  which  adaptation  is  accomplished. 
Items  of  interest  in  this  program  were  the  time  required  to  adapt  the  air- 
frame from  its  basic  un-augmented  condition  to  damping  ratio  =  0.7,  and  the 
steady  state  value  of  damper  gain  upon  completion  of  adaptation.   These  values 
are  listed  in  Table  XI  for  each  flight  condition.   Other  items  listed  are  the 
theoretically  predicted  gain  for  adaptation  and  the  basic  airframe  short  per- 
iod damping  ratio.   Fgiure  30  shows  that  the  steady  state  damper  gain  exhibit- 
ed small  amplitude  fluctuations  about  a  mean  value  during  operation,  and  it  is 
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Figure  30.   Time  History  of  Airplane  &  Control  Surface  Motion 
and  Pitch  Damper  Gain  During  Adaptation  Process  With  Adapt- 
ation Servo  Lock 

Initial  K^  =  0,   Adaptation  Servo  Gain  =  1. 

Flight  Condition  7  -  M  =  1.8  ©  25,000  ft. 
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Table  X 
of  Freedom 
Servo  Lock 
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Listing  of  Computer  Output  from  a  Three  Degree 
Flight  Simulation  With  Adaptation  Loop  Using 
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Table  XI  -  Comparison  of  Steady  State  Damper  Gain  Generated  by  the  Adapta- 
tion System,  with  Theoretical  Gain  Required  for   =  0,7  at 

ch  Flight  Condition,  tin?  to   edapt,  &  Basic  Airframe  Damping 
Ratio  also  Tabulated. 


Adaptation  Servo  Gain  1.0 


Flight 
Condition 


Rasic 
f 


Time  to 
apt 


Theoretical 

Kd 


Experimental 
Kd 


1   H  =0.9 
10,000  ft 

0.440 

0.156 

0.16 

2    M  =  ! .  2 

10,000  ft 

0.177 

2.2 

0.303 

0.29 

M  *  1 . 
10,000  ft 

0.156 

2.1 

0.284 

0.28 

M  =  0.9 
25,CCC  ft 

0,244 

1.7 

0.249 

0.26 

5   M  -  1.2 

25,000  ft 

.  03 

2.5 

0.372 

0.36 

6   M  ■  1.4 
25,000  ft 

0.105 

2.5 

0.384 

0.38 

7   H  v  1.8 

25,000  ft 

0.096 

2.5 

C. ' 

0.34 

8   M  =0.9 
45.000  ft 

0.1  ' 

2 .  2 
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0.  \\ 

9   M  =  1.2 

45,000  ft 

0 .  066 

o.c 

0.56 

10  M  =  1.4 

45,000  ft 
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2.9 

0.503 

0.51 

11  M  -  1.8 
45.000  ft 

0.058 

■  .„_ 

2.8 

0.467 

0.47 
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It  is  seen  that  the     tation  times  compare  favora';      th  those  of  the 
MIT  system,  and  as  would  be  expected,  those  flight  conditions  with  lower 
values  of  basic  damping  ratio  exhibited  longer  adaptation  times.   The  ex- 
perimental values  of  steady  state  damper  gain  could  be  determined  at  best 
with  two  place  accuracy,  but  even  so,  reasonable  agreement  with  the  theo- 
retical values  was  obtained. 

Considerations  in  the  Practical  Mechanization  of  an 
Airborne  Digital  Damping  Ratio  Computer 

Perhaps  the  best  method  of  reviewing  the  finally  "optimized"  Adapt- 
ation process  is  to  follow  through  the  flow  diagram  of  Figure  31.  which 
illustrates  the  program  sequence  to  be  followed  by  the  proposed  airborne 
damping  ratio  computer.  The  primary  difference  between  the  airborne  com- 
puter and  the  §-uJn  solution  programmed  into  the  three  degree  of  freedom 
flight  simulation  is  the  requirement  for  an  analog  to  digital  converter  for 
pitch  rate  input  signals  and  a  digital  to  analog  converter  for  output  signals 
to  the  adaptation  servo.   It  was  found  from  the  simulation  program  that  simple 

trapezoidal  differentiation  of  0  and  6  yielded  sufficiently  accurate  values 

• »     ...  . 

of  9  and  G  for  use  in  the  £-%  calculations. 

From  a  short  preliminary  survey  of  computer  requirements  it  was  deter- 
mined that  a  special  purpose  digital  computer  having  a  130  word  storage 
capacity  and  a  repertoire  of  12  arithmetic  and  manipulating  instructions 
would  suffice  for  this  application.  And,  if  it  is  assumed  that  each  in- 
struction requires  5  microseconds  for  execution  (a  representative  state- 
of-the-art-estimation),  the  sequence  from  raw  data  input  to  generation  of 
a  new  adaptation-servo  command  will  require  12,600  microseconds,  or  in 
other  words,  the  total  cycle  will  have  a  maximum  possible  repetition  rate 
of  83  times  per  second.   However,  the  cycle  utilized  in  this  study  was  re- 
peated only  40  times/second,  which  leaves  a  comfortable  execution-speed 

margin. 
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Figure  31.   Flow  Diagram  of  the  Proposed  Damping 
Ratio  Computer  for  Use  in  the  Adaptation  Loop 
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VIII  -  CONCLUSIONS 

It  has  been  demonstrated  that  a  simple,  special-purpose  digital  com- 
puter can  be  utilized  to  determine  the  natural  frequency  and  damping  ratio 
of  the  short-period  mode  of  the  free  longitudinal  oscillations  exhibited 
by  a  representative  high  performance  fighter  aircraft.   And8  it  has  also 
been  shown  that  the  calculated  values  of  damping  ratio  can  be  used  in  a 
closed  loop  manner  to  vary  the  gain  of  a  self  adaptive  pitch  damper  which 
is  capable  of  maintaining  constant  airframe  short-period  damping  ratio  over 
the  full  speed  and  altitude  range  of  the  vehicle.  The  system  investigated 
is  capable  of  performing  the  adaptation  process,  from  basic  unaugmented 
airplane  to  £  =  0.7,  in  a  period  of  time  ranging  from  1„4  to  3.1  seconds 
depending  on  flight  condition;  and  this  speed  of  operation  was  shown  to  be 
similar  to  that  of  the  most  recent  MIT  flight-vehicle  self-adaptive  system. 

In  conducting  an  evaluation  program  with  the  proposed  mechanization, 
in  conjunction  with  a  longitudinal  three  degree  of  freedom  flight  simula- 
tion, it  was  found  that  measurements  of  pitch  rate  rather  than  pitch  angle 
were  necessary  to  extract  the  short  period  characteristics  from  the  combin- 
ed phugoid  and  short  period  motions.   First  and  second  derivatives  of  pitch 
rate  are  then  determined  in  the  adaptation  computer  and  sampled  values  of 
8,  9,  and  6  set  into  the  iterative  routine  for  determination  of  short  period 
damping  ratio  and  natural  frequency.   The  system  requires  no  perturbation 
signals  or  limit  cycling  for  operation  of  the  identification  process,  how- 
ever, the  system  will  adapt  only  while  the  airframe  exhibits  pitching  motion. 
It  should  be  emphasized  that  the  mechanization  investigated  in  this  study 
is  not  a  complete  longitudinal  control  system,  in  that  there  is  no  pro- 
vision for  external  pilot-input  signals.   However,  it  was  proposed  in 
Chapter  V  that  Corbin's  forward- loop -gain  computing  system  could  be  con- 
verted from  analog  to  digital  form,  and  used  in  conjunction  with  the  present 
computer  routine  to  enable  the  introduction  of  external  control  inputs. 

89 


tg  . 

study  of  of  forcing 

the  convergence  of  t,~Mn   iterative  c.        is.   It  is  felt 

that,  if  the  £-U)n  calculations  could  be  n&de   100  percent  convergent^  the 
speed  and  accuracy  cf  the  adaptation  process  could  be  greatly  enhanced. 
Other  factors  which  would  require  investigation  before  a  practical  mechani- 
sation could  be  accomplished  are:  1.   the  effect  of  noise  (atmospheric  tur- 
bulence) introduced  into  the  forward  loop,  2.   the  effects  of  control  sur- 
face rate  snd  amplitude  limiting^  and  3.   the  influence  of  transducer  in- 
accuracies and  noise  on  the  operation  of  the  adaptation  loop. 
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Appendix.  I 

Table     of  Stability  Derivative  5   of  a, 
f\e pre-sentai'ive     Hian     Herformtxnce. 

Fiojn  ter  "At  totcK     Aircraft 

Altitude        10,000  ft. 
M^O.B  M~  /.  2-  M-J.4-  M-/3 


Om;± 

—,<s~ra 

— .  4-4-/ 

-.333 

Cm  $ 

~/o.B3 

-4-.3S& 

-3.2.+<h 

C^c, 

-.  +  7S 

-.  G6S 

-  .S'/S 

C"cz 

-.  83/ 

+  2..03 

+-  .73S- 

Cl« 

4-.Z.2. 

<h.J8 

£..  8<?S- 

G-* 

.S03 

.  2-73 

.  2.14- 

Cn# 

.ZOS 

.138 

.112. 

-.34-4- 

-293 

-  .22-4- 

C"<fcc 

.0/34- 

.  0031 

.002.  9 

0 

Onr 

-.4-31 

~.+4-4- 

-.  3  GO 

Onp 

.oao 

.  O  6  4- 

.01T4- 

CjR 

-.130 

-.  /2.9 

-.117 

C*6» 

.HI 

.OS3 

.08S 

C*Ju 

.  00  a 

.OOGS 

.0037 

CjLr 

./£/ 

.  1+4- 

.  113 

• 

CjLp 

~.3+S 

-.2.SS 

-.2.1 

*  * 

&ye 

-J./S 

-  /.Of 

-  .  9/             -.  .. 

c* 

-' 

cv^ 

.  720 

.S73 

.  43S 

Cy» 

cD. 

* 

.  J 

All  derivatives    per    radictn 
or    rocdian per*  ^second 
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ia.b\e    of    utccbi/iiy    Derivatives 


Altitude      2.S,ooo    ft. 


M'O.Q 

M*  J.  2. 

M~  /■*- 

M~  /.a 
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-S3J 
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.3SO 

■    T 
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.     * 

.  jr          1 
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I cxble    of    otoLbil'ity    LJerivaUves 


Altitude       ?$tOOO    ft. 


m-  o.e 

M~  J.Z. 

M^  /•*- 

M~  /•£ 

CM}t 

~.99J 

-.  7+5 
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-  .3€7 

(sfilq 

-4:33 

-4-.  7/ 

-3.83 

■-  a.eo 
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-  .&02. 
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-  J.3G 

•h  2..3T 
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4-.74> 
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3.4-3 

2.4-7 

c*,. 

.  72.5 
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.352. 
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.  267 
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.zoe 
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^n6r 

-.335 

-.3S4- 
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-.2.0/ 
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Cr)r 
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./OB. 
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/  34- 

./£S 
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.087 

Cj.Su 
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.0/2.7 

.0038 
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CjLr 

.  2.4-2. 

.J7& 

.13d 

.//4- 

Cj.p 
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-.2.6^ 

Cye 

-A  30 

-/./  3 

-  J./O 

-.57       ' 

Cyr 

-■  ^. 

C^Sr 
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•5€  8 

.4-13 

CyP 
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APFEl 
Table  of  Airj  al  Characterise 


Airplane  mass 


xz 

Wing  area 

i-Jing  span 

Wing  mean  aerodynamic  chord 

Wing  aspect  ratio 

Xe 


1, 110  slugs 

.42   slug-ft. 
129,    533  slug   f 
142,398   slug-ft62 
8,219  slug-ft/ 

394.84  ft.2 

36.58  ft, 

12.24  ft. 

3.39 

13,95  slug-ft. 

1  s.259  rad/secs  (Cruising) 
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A  FM  .  I7ES 

DERIVATIVE  CO; 

C...  Rolling  moment/  aileron  deflection  -  little  imports 

in  dynamic  stability. 

C  Rolling  moment/  rudder  '        n  -  cross  coupling  te 

of  little  importance  in  dynamic  stability, 

C,  Rolling  moment/  roll  rate  -  viscous  damping  in  roll, 

balances  aileron  moment  to  produce  a  constant  rolling 
velocity  for  a  given  aileron  deflection;  has  small 
damping  effect  on  lateral-directional  spiral  mo 

C.  Boiling  moment/  yav?ing  velocity  -  contributes  slightly 

r  to  Dutch  roll,  most  important  cause  of  spiral  divergence 

mode. 

C.  Rollimg  moment/  sideslip  angle  -  dehedral  effect,  large 

*  negative  value  will  produce  good  spiral  characteristics 

and  poor  Dutch  roll  characteristics. 

C-  •  Rolling  moment/  sideslip  angular  rates  usually  negligible 

in  magnitude. 

C  Pitching  moment/  elevator  deflection,  control  derivative, 

i  little  importance  in  dynamic  stability. 

C  Pitching  moment/  pitch  rate,  viscous  damping  in  pitch, 

q  large  effect  en  damping  of  longitudinal  short  period 

mode. 

C  ,  Pitching  moment/  rate  of  change  of  angle  of  attack^  due 

"1<3:  to  aerodynamic  time  lag  and  aeroelastic  effects^,  contributes 

to  damping  of  longitudinal  short  period  mode. 

C  Pitching  moment/  angle  of  attack  -  produces  restoring  moment, 


v.: 


a  establishes  frequency  of  longitudinal  short  period  mode. 


C  Pitching  moment/  sideslip,  angle  -  aerodynamic  cross  coupling 

P  derivative,  usually  r     ible  in  magnitude. 

C  Yawing  moment/  rudder  deflection  -  control  derivative,  litt ' 

"<*  r  importance  In  dynamic  stability. 

C  Yawing  moment/  yaw  rate  -  viscous  damping  in  yaw,  major  contri- 

r  butor  to      ng  in  Dutch  roll  (lateral-directional  short  peri- 
od mode) . 


Ill 


C  -  aerodynamic  cross  coup  I1 

p  ites  slightly  to  damping  in  Dutch  roll* 

C  Ya«.        it/  sides!     ;le  -  produces  directional  re- 
^  nt,  establishes  frequency  of  Dutch  roll  mode. 

C  .  Yawing  moment/  sideslip  angular  rate  -  contributes  slightly 

"&  to  damping  of  Dutch  roll. 

C  Yawing  moment/  aileron  deflection  -  control  cross  coupling 

<^i  derivative,  little  importance  in  dynamic  stability. 

C  Side  force/  sideslip  angle  -  contributes  to  damping  of  Dutch 

fi  von. 

Cy  Side  force/  roll  rate  -  aerodynamic  cross  coupling  deriva- 

~p  tive,         j  usually  negligible. 

Cv  ie  force/  yaw  rate  -  magnitude  usually  negligible. 

r 

C  .  Side  force/  sideslip  angular  rate  -  magnitude  usually 

0  negligible. 

C  Lift/  angle  of  attack  -  effects  frequency  of  longitudinal 

JC£  phug  i   ode,  has  major  effects  on  aircraft  maneuverability. 

CT  .  Lift/  rate  of  change  of  angle  of  attack  -  due  to  aer-odynamic 

"GJ  time  lag  and  aeroelastic  effects,  has  negligible  effect  on 
dynamics. 

C  Lift/  elevator  deflection  -  magnitude  negligible  for  conven- 
tional aircraft,  has  effect  on  tailless  aircraft  and  those 


\ 


with  short  tail  length. 


C  Airplane  drag  coefficient  -  primary  contributor  to  damping  of 

longitudinal      Ld  osc      on, 

C.  Drag/  angle  of  attack  -  negligible  effect  on longitudinal 

00  phugoid  mode  because  this  mode  takes  place  at  constant  angle 

of  attack. 
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TX  IV 
EQUAT  i        [OTION 
The  equatic  -  derived  from  tho  t: 

tons:   Force  -  Mass  x  Acceleration  a. 

Moment  =  Moment  of  Inertia  x  Angular  Accelerate 
From  classical  mechanics,  the  equations  of  motion  of  a  rigid  body  in  free 
space  with  no  aerodynamic  or  gravity  forces  acting 

f~X    =   rn(u  +  qw  -rv) 

Fy  -  m(  v  +  ru   -pw) 

F%  -  m  ( iv  + pv  -  <£  u ) 

L=     ph  -K}r(h-Iy)  "  IteCp^+ir)  -hIxy(rp-$)  +  Iy£(r*-<f) 

M^   yly  +  (p%-r)Iyi-(p+yr)Ixy  +  rp  (Jx~Il)+  (pz-r*)IXi 

N=  rli    t  (^r-/>)7Xi  -(y+rp)IYi  +  p^(Iy-Ix)  +  (<f-p*>IXy 

Note  that  all   symbols  are  defined  in  figure   1.      It   is   seen  tl  e   last 

t^o  terms    in  each  of  the   first   three  equations  are  thi  of  an  angu- 

lar and  a   linear  velocity  and  thus  represent  centripet.  t rat  ions  dir- 

ected along  the  re  we  axes.      The   last    four  che   last 

three  equations   represent  moments   created  by   inertial   cross   coupling,,      These 
terms  arise   from  a  vai  s    including  g  1   torques 


and  i 

ertial  cc      ' 

pointed  out  that  times  i  trous 

results,  when  the  arnplitu        rates  of  aircraft  motion  be<  [ref. 

3).   It  is  iircraft  posses  lich 

causes  the  products  of  Inertia  I       I   to 

inertial  coupling  tei  oro  each  e<; 

force  due  to  gravity,  f< 


P       »  , 

rnu    =     rv  -ow  -  met  sin  9  -  Cp-p-V^S   1~  Engine  Thrust 


9       L 
mv  -  pw  - ru    +■  ma  cosS  s'm0   +  Cy-£  V   S 

mvs/  =  qu    -pv    -mo    cosO  cosy)    'CiT'V^S 


$Iy  =  f»(Ii-Ix)  t-r^xz  -P^hl  -heaJer  +Cm-j-V^Sc 


P       ? 
rl^p^dx-Jy)  +pIXz    -<j.r  IXii'Ixeu}e%-i'Cr>T  V   ^ 


The    E      Cl>        terms    represent  parts, 

Xe 

which  when  multiplied  by  a  yaw  oi  xt   velocl          s  a  gyro- 
scopic processional  torque  about  the  pitch  aw  axis  re            These 
terms  are  significant  in  a  turbojet  powered  aircraft.   It  should  be  noted 
that  the  wingspan,  b,  is  ients  in 
the  roll  and  yaw  equations  whereas  mean  ".          used 


mome] 

A 
comi    ■ 
into  rnc  ies, 

1/  »     V   ;         airo/cLne      forward      velocity 
\/v    V'6in£       «    l/# 
W&  V sin  (Z    *•    V 'a. 

v  «  yi 

w«.   Vex 

Whera  cc    1.^5  ,  the  aerodynamic  at  ;gless 

velocity  vector  to  the        -axis.   The 
complete  equations  of  motion  with  the  aerodynam<        -tents  it      ded 
form  are  shown  in  table  T,  The  reader  will  see  t       aerodynamic  rc- 
tary  derivative  terms  are  multiplied  by  certain  add      1  flight  condit 
parameters  in   ©rde,        e  value  of  the  derivativ 
as  possible  over  a  wide  range  of  flight  conditions .     example  is 
of  the  term  -£-  wich  derivatives  taken  with  res* 
terms  for  other  derivatives  are  listed  in  thts  table  of  symbols. 

One  final  set  of  equations  must  be  used  to  ei  on  of 

angular  rates  abc*  s  to  angular  rates  with  res        he 

earth  coordinate  s-       The  necessary  Euler  equations  are       5.)". 


115 


Q  =     q  co  3  cf>    —  r   ■sin  <£> 

<fi  =     p    +  <q   Glnd)    torn  <p  +  r  coj<p    tan  S 

9  -      (%   s)n<t>   +  r  cos  $ )  s<so  & 

•  m 

G  and  %     can  be  integrated  to  obtain  0  and  0  which  are  used  in  the  equa- 
tions of  motion. 
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V 
DERIVATION  OF  THE  LO:        tt,  SHORT-PERIOD  TRANSFER  FUNCTIC        HE 

EQUATIONS  OF  MOTK 
Of  the  two  longitudinal  modes  of  airplane  motion,  the  phugoid  and 
the  short-period,  only  the  latter  mode  Is  of  interest  in  stabil        cita- 
tion considerations  because  the  phugoid  normally  does  not  cause  any  deter- 
ioration of  handling  qualities.   In  Section  III,  a  transfer  function  of  the 
mode  of  interest,  in  stability  augmentation  systems  must  be  included  in  the 
block  diagram  of  the  proposed  fitch  damper  to  enable  an  analysis  of  the 
system.  Thus,  the  longitudinal  short-period  transfer  function  must  be  de- 
rived from  the  airplane  equations  of  motion,   The  legrees  of  freedom  of 
interest  being  pitch  and  normal  translation  for  which  the  corresponding 
equations  are  shown  below, 

25     AC      ^  -  Ac2-    ^  A.      Ac     r         a„ 

6   =  -j^—  Cm',t    it   +•  Zy?y   "2^  Q        1y  Cm<*  AC° 

CT,=G+    ^     COSO    -^rC^CC 

The   pitch  damping  derivatives,  Cm  C/tl^        are   usually  added  together 

to   form  an  equivalent  Cmo^o,       >    and  engine   torque  and  tail    lift   terms 
have  been  neglected  because  oC  their   low  mangitudes.      Ass  neces- 

sary  for  the  derivation  of  the  transfer   function  are:      1„  illibriu 

condition  exists  al  he  2     ixis  when  airplane  weight  =  ait  lift, 

W=*   4-SVzC      Ortnm  I.    cos  6^1. 

The  difference  between    CC    an  <X  -trim  will   be   denoted   Cc, 


11? 


?he   equations   of  motion   then  bee 


ome 


& =       ly     ^    '*    *"     £Wy~  C^?^  ^   *     77"  C/wo,   <^/ 


A 

OL=    9    ~     -JffiT    CLnr.   ^i 


The  next  step  is  the  introduction  of  the  differential  operator^  S,  and  the 

Q 

elimination  of  angle  of  attack  from  the  equations  in  order  that  the  — I — 
equation  may  be  obtained. 


SL&  -  i^Cm)t>t  +  j^r  Cm^k  66  -h  j—  Cnn^  00, 


ocs  =es  -  j^j  cLql  cc, 


or:  CC,   =  9S 


^  *"   TZU'Cl 


mv  ^Lcc 


Substituting     CC, 


cl       Ao2-^  co_    Ac    ^  69 _       Ac 


ff)V  ^Loo 


Q 

The    -s —     transfer    function   is 

't 


Ac  /  fy  * 


>t    "      c/sL+c-A-r    -A£±r       r        AC2-  r         A     r  AG  r         \ 


Rearranging 


e 


it     s(^+skkcLa-ggc„f)-(A%cm^cLa  -4&-C„J 
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The 

9    _       fo(Si-Te) 

re: 

Kq  '  vaness 

Te  Light  path  fii.ie  constant 

£      the  short-period  d i      ratio 
&//,   Ls  the  short-period  natural  frequency 
It  can  be  seen  that  the  steady  state  transfer  function  is  essentially 
a  pitch  rate  command  system  for  a  given  horizontal  s     Izer  deflection. 
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IX  VI 

[C  SYSTEM 
Although  analog  computers  are  generally  used  to  simulate  system 
dynamics,  the  systems  Investigated  in  this  study  were  simulated  on  a  stored 
program,  genera   •  urpose  digital  computer  using  techniques  similar  to  ana- 
log computer  programs,   The  techniques  are  similar  in  that  both  methods  use 
some  form  of  integration  combined  v;ith  other  operations  such  as  summation, 
subtraction,  ".[.on  to  solve  a  set  of  simultaneous  different 

equations  of  motion;  each  equation  of  motion  representing  a  degree  of  free- 
dom of  the  system.   An  analog  conputer  program  may  be  regarded  as  a  cont * 
ous  parallel  solution  of  the  equations  of  motion,  whereas  the  digital  method 
employs  a  step  by  step  serial  solution.   A  number  of  digital  integration 
techniques  are  available,  ranging  from  simple  rectangular  integration  to 
highly  accurate,  but  time  consuming  methods,  such  as  the  Runge-Kutta  numeri- 
cal integration.   The  technique  selected  for  this  study  is  a  modified  form 
oi:  the  trapezoidal  rule  and  is  shown  below. 

•'here:  n-k  denotes  values  computed  in  previous  cycles 

h  is  the  width,  in  time,  of  the  computation  interval, 

A  major  difference  between  the  digital  solution  and  the  concept  of  the 
analog  Clow  diagrai  1 ies  in  the  fact  that  the  digital  method  performs  in- 
tegration of  lower  ordered  derivatives  first.   As  an  example,  the  steps 
in  the  computation  cycle  of  a  second  order  system  solution  are  shown  below: 

%n  '   X«-/  +■-]*-£  3Xn-,  "  X/rzJ 

X0*-d$WnX„~u;fX„ 
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The  computation  cycle  Is  repeated  it   high  rate  and  the  values  of  X  X,  £  X 

•rated  in  each  pass  re]  uts  on  a  time  history  of.  system  motion. 

TJhen  a  more  i      :•:  system  having  «  greater  number  of7  aquations  of  motions 

Ls  simulated,  the  integrations  for  e."»ch  equation  of  motion  are  performed 
in  a  serial  manner  such  that  one  pass  is  executed  through  all  equations 
before  recycling. 

The  accuracy  of  this  technique  is,  of  course,  a  function  of  the 
vidth  of  time  increment  between  computing  cycles,  and  it  was  noted  that 
when  using  the  method  described  above,  integrator  errors  were  manifested 
as  sv,   undamping  of  the  system  dynamics.   In  some  lightly  damped  cases  it 
was  found  that  cumulative  integrator  errors  actually  caused  the  solution  to 
diverge,  the  corrective  action  being  the  use  of  finer  calculation  inter- 
vals.  It  has  also  been  noted  that  the  required  fineness  of  calculation 
intervals  is  a  function  of  both  the  period  of  the  system  mode  of  highest 
frequency  and  the  damping  of  this  mode.   A  rough  rule  of  thumb  is  the  use 
of  200  divisions  per  period  for  £   =  0.1,  down  to  50  divisions  per  period  for 

<7  =  .5,  to  obtain  one  percent  solution  accuracy. 

The  following  five  pages  are  sample  listings  of  two  of  the  FORTRAN- 
programming- language  programs  written  for  this  study.   The  first  is  a 
second  order  system  and  the  second  a  three  degree  of  freedom  flight  simula- 
tion.  The  programs  were  run  on  both  the  CDC  1604  and  the  IBM  7090  general 
purpose  digital  computers. 
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SECOND  ORDER  PLANT  IDENTIFICATION  ROUTINE  3 

_1^£AD  IMEUJ    TAPE  3,2,  ACHM...ALT,  DENSTY,  Vt:LX,  DYNPRS,  THETA.  CTMET 
IAjSTHETA,  ALPHAT,  THRUST 

2  FORMAT     (FA. I,     16,     F3.7,     F7.2,     F8.1    ,    4F7.6,     F8.2) 

KLAD    INPUT    fAPB    3,3,    CMIT,    CMUA,    CMA,     CLA,    CLUELE,     CUU,     TIIEDD 

3  FORMAT  (6Fo.3,  P8.4) 

WRITE  OUTPUT  TAPE  2, A,  ACHM,  ALT 

'"^m^-I"  -4^/58H  --r-HKEE  DEGREE  ..OF..  FREEDOM.  FLIGHT  SIMULATION   MACH 

lNUMBER*  F4.1,  12H   ALTITUDE=  16) 
IN  I TXAL- CONDITIONS 

J  =  i 
J  1=1    

K  =  l 


50 


Ll=l 

TIME=0. 

SRVGN=.25 

-DAMPR=.7 

DMPGN=0. 

-DAMPR2  =  .4 

ACCLN=0. 

OMEGAN=0. 

ALPHA=.5 

ALPHD=0. 

ALPHD.=0. 
TH£TD=G. 

THETDi=0. 
-  -THEDDi  =  THEDD 

THTTD=0. 

THT.TD2  =  0. 

.VELCTY  =  VELX 
— vfi  vn=n,  > 


VELXDi=0. 
_ACTPOS=0. 
ACTVL=0. 
ACTVL1=0T 

ACTCL=C. 
-ACTCL1=Q.._ 
SIGNAL=0. 


ALPHA2=.5  

ALPHD2=0.  .  „  - 

__  THETA2  =  .5  

THETD2=0, 

THEDDl^THEDD   

5  WRITE  OUTPUT  TAPE  2,6 

6  FORMAT  {105H    TIME-SEC       THETA      THETADOT    FWD 
1UAT0R  POS  NORM  ACCEL   DAMPER  GAIN  OMEGAN   '     DAMPR) 

7  WRITE  OUTPUT  TAPE  2, 8, TIME,  THETA, 


1DMPGN,  OMEGAN,  DAMPR 

8  FORMAT  (9F12.'6) 

IFU_-1000)9, 12*12 

9  L=L+I 

IF(Ll-26b)51, 52,52 
1  L1=L1+1 

GO  TO  14 


VELOCITY  ACT 
MPR) 
THETD,VELCTY,  ACTPOS,  ACCLN, 
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SECOND  ORDER  PLANT  I  DENT  I  PI CA1 ION  ROUTINE 


52 
12 

14 


15 

16 

21 

27 
30 

31 

25 


LI 

GO 
CA 
/\C 
AC 
AC 
AC 
AC 
AC 
IN 
VE 
VE 
VE 

15  9 
Vc 

.CA 
DY 

.  AC 
AC 

_AM 
AO 
SO 
TH 

_IH 
TH 

__TH 
TH 

1PH 
AL 
AL 

v  AL 

.  CA 
AC 

.  PI 
SI 

j:TI 
DI 

_JH 
AD 
DM 
IF 
J  = 
GO 
J  = 
CA 
DA 
GA 
BE 
OM 
DA 
IF 
Jl 
GO 
Jl 
IF 
IF 


_r_a  io 

LL    EXI 

TUATOR 

ppqs«a 

TVL1=A 

TVL=AC 

TCLI=A 

TCL=-A 

TEGRAT 

LX=VEL 

LX01=V 

LXD=(T 

JL'i.2  9) 

LCTY=V 

LCULAT 

N.PRS  =  D 

IYY=DY 

2VIY=( 

V  =  DYfviP 

V£RM=D 

LUTION 

GTA=TH 

ETD1=T 

ETD=TH 

EDDI=T 

EDD=(  A 

AT)  ) 

PHA=AL 

PHDl.=  A 

PHD=TH 

LCULAT 

CLN=(C 

TC.H    DA 

GNAL=T 

ME=TIM 

FFEREN 

TTD=(T 

APTA7 I 

PGN=DM 

( J-5)l 

J+l 

TU    14 
1 

LCULAT 
MPRlrD 
MMA=DA 
TA-SQR 
EGAN=- 
MPR1=- 
(  J  1-2  0 
=J1+1 

TO  2i 
=  1 

(OMEGA 
(OMEGA 


T 

SIMU 
CTPOS 
CTVL 
TVL+. 
CTCL- 
CTVL* 
ION  0 
X+.0  0 
ELXD 
HRUST 
)-(32 
ELX 
ION  0 
ENSTY 
NPRS* 
ACIYY 
RS/(  V 
YNPRS 

OF  E 
ETA+. 
HETD 
ETD+. 
HEDD 
CI  YY« 


LATION 

+.00  23*(3.»ACTVL-ACTVL1) 

0025*(3.*ACTCL-ACTCL1 ) 


49.-ACTP0 
F . LONITUD 
25*(3.*VE 

/  1 1 1  0  .  )  -  { 

._18*S1NF( 

F  AERODYN 
*( VELCTY* 
12.24/i29 
*12.24)/( 
ELCTY*111 
/  1 1 1  0 . 
QUATIUNS 
0025*(3.* 

002  5* ( 3. * 

CMIT*ACTP 


S*l 

INA 
LXD 

AOV 
THE 


225.+SIGNAL*1225. 

L  ACCELERATION  TO  OBTAIN  VELOCITY 
-VELXD1 ) 

I KM*CDO)-(AQVERM*(  (CL A* ALPHA) **2.  )/ 
T  A  )  )    


(3.14. 


C  PARAMETERS  DEPtNDENT  ON  VELOCITY 
)*394.34 


AMI 

*2. 

532. 

2.*VELCTY) 

0.  ) 


OF 
THE 


MOTION 
TD-THETD1) 


THEDD-T 
US) 


+  UC2VI  Y«CMQA*THETD) +  (  ACIYY  *C  MA*  (  ALPHA-AL 
HD-ALPHDi ) 


HETA)/VELCTY)-(AMV*CLA* ALPHA) 

T  I  0  N    L^__L 

(1110.*32.18) 


PHA+.0025*(3.*ALP 

LPHD 

ETD+( 32.18*C0SF(T 

E  NORMAL  ACCELERA 
LA*ALPHA*DYNPRS) / 

MPER  SIMULATION 

HETD*DMPGN 

E+.00  5 

TIATE  PITCHING  ACCELERATION 

HhDD-THEDDl)/.005 

ON    SERVO 

PGN+(  (  .7-DAMPR)  *  SR\/GN*  .005  ) 

5,16,16 


E  NATURAL  FREQUENCY  AND  DAMPING  RATIO-  ASSUMED  DAMPR=.4 

AMPR2 

MPR1*THE602/THETD2 

TF(GAMMA**2.-THTTD2/THETD2) 

GAMMA+8ETA 

(THT'TD/(2.  *OMEGAN*THE_DDJ  )-(  (  GMEGAN*THETD  )  /  (  2  .  *  THE  DO  )  ) 

)20,31,3i 


N-8. )25,26,26 
N)  26, 32, 32 
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SECOND  ORDER  PLANT  IDENTIFICATION  ROUTINE  3 


32  IF(DAMP;<1-1.  )33,26,26 

33-1  FID AMP.Kl 126, *34,  34        

34  DAMPR=DAMPR!  "- 
DAMPR2=OAMP_Rl 
GO  TO  29 
26  DAMPR1=DAKPR2 

3  5  GAMMA=DAMPRi#THEDD2,/THETD2 

BET A=  S w RI F  { G A M MA * *2  .  - T H T T D 2 / . T H ET D 2  ) L 

OMEGAN=-GAMMA-B£TA 

DAMPRl=-(THTJD/(2.*OMEGAN*THEDD)  )-(  {  OMEGAN*THETD  )  /  (  2.  *THEDD  )  ) 

36  IF(Ji-20)37,38,38 

37  Jl=Ji+l 

GO  TO  35   N  '; 

3B..jiJLi  "  : . 

I  F  (O.MEG AN- 8.  )  3  9,40,40 
39  IF<0MLG'AN)40,nl,4I 

41  LF(DAMPR1-1. )42,40,40 

42  IF (DAMPR1 ) 40, 34, 34  .  ^ 

4  0  DAMPR=.7 

■ WRITE..  OUTPUT..  TAPE  2,43  

'43  FORMAT  (31H   NO  SOLUTION  FOR  DAMPING  RATIO) 
29  THETD2=THETD 

THEDD2=THEDD 

THTTD2=THTrD 

IF (K-iO) 17,18,18 

1 7  .  K  =  K  +  1   : 

GO  TO  7 

18  K  =  l 

GO  TO  5 


E  \  D  (  I ,  J  ,  0  ,  0 ,  0  ,  0 ,  0 ,  0 ,  C  ,  0 , 0  ,  0 ,  0 ,  0 ,  0  ) 
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PLANT  IDENTIFICATION  ROUTINE  SECOND  ORCER  2 

XEC      ___ 
1  =  1 

J=1    

K=1 

XBO=0« 

XDD0=C. 

XC  =  1.  l 

XDN=Q. 

XN  =  1.  . 

TIME=0« 

XDCO=C.         

XDDOC=Q. 

WRITE    OUTPUT    TAPE    2,1  

2    cT/pRl"H1//kkH    PLANT    IDENTIFIC^^^OuflNT-^SECOND    ORDER) 

XDDN=-2.i|«XDN-36o*XN 

XDINCR={XDDN+XCD0)*.C0D25 

XCN=XCO+XDINCR 

XCDO=XCDN 


XINCR= (XDO+XDN )*. 00^25 

XN=XN+XINCR 

XDO=XDN 

TIME=TIME+.0035 

IF(I-1003)U,5,5 

1=1+1  

GO  TO  2 
1  =  1 


WRITE  OUTPUT  TAPE  2,6  ~ ' 

'Si!    MH    //55H 11HE ^CtL£RA.UON^VELOCITY„         DISPLA 

^I.TE_CUTPI^I  TAPE  2*7»  TIME,  xnnu,  xdn.  xm 


7  FORMAT  (UF13.6) 

IIE^lDTTK,rAh*IJLATI0N  °F  NATURAL  FREQUENCjr_A_ND  DAMPING  RATIO  GIVE 

A  STARTING  VALUE  OF  DAMPING  RATIO  =  rt.ii      •■ 

WRITE  OUTPUT  TAPE  2\8  

8  FORMAT  (25H      CMEGA-N       DAMPRAT) 

9  ALPHA=DAMPR#XDCO/XO 
BETA=SGRTF(ALPHA«*2.-XDD00/X0~j"  ~*-"~ 

ovegan=-alpha+beta        

omeg^bIf^ 

IFIOMEG-100. )  15,16,16  ~^~^ 

15  DMP=A8SF(DAMPR)  m 

IFIDMP-1.  )17,16,16 

16  WRITE  OUTPUT  TAPE  2,18 ■ 

18  FORMAT  <19H  SOLUTION  DIVERGING) 

DAMPRs.4  .i' 

J  =  1 

20  ALPHA=CAMPR«XD0O/X0  

BETA  =  SCRTF(A'LPHA**2.-XCDC0/X0)  ~ 
OPEGAN=-ALPHA-BETA  

U"  tb-Adi>h(0  MEGAN) 

IF(0MEG-1C0. )21, 22,^2  ' 

21  DMP=ABSF(DAMPR)    

IF(DMP-1.)23,22,22  ~^T 


PLANT  IDENTIFICATION  ROUTINE  ..SECOND  ORCER  '2 
22  WRITE  OUTPUT  TAPE  2,24 


2U  FORMAT  (19H  SOLUTION  DIVERGING) 
GO  TO  12      : 


23  WRITE  OUTPUT  TAPE  2,25,  OMEGAN,  DAMPR 

25  FORMAT  (2F13.6)  

IF(J-20)26,12,12 

26  J=J+T  _^ 

GO  TO  20 

17  WRITE  OUTPUT  J APE  JZs .10  , ..J0MEGAN,__O.A.MPR. 

10  FORMAT  (2F13;6) 

I FU-2C)  11,12,12 

11  J=J+1 

GO  TO  9  = 

12  J  =  l 

xo=x'N       : 

XDOO=XCN 

X€DOO=XDDN    

IF,(K-1»G>13,14,14 

13  K=K+1  

GO  TO  2 

1U  CALL  EXIT  


END(1, 0,0,0,0,  0,0, 0,0, 0,0,  CO,  0,3) 
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^Lor  investigation  of  a  ^adap 


A  simulator 
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